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ABSTRACT 

Aims. In this work, we aim at providing a consistent analysis of the dust properties from metal-poor to metal-rich environments by 
linking them to fundamental galactic parameters. 

Methods. We consider two samples of galaxies: the Dwarf Galaxy Survey (DGS) and the Key Insights on Nearby Galaxies: a Far- 
Infrared Survey with Herschel (KINGFISH), totalling 109 galaxies, spanning almost 2 dex in metallicity. We collect infrared (IR) to 
submillimetre (submm) data for both samples and present the complete data set for the DGS sample. We model the observed spectral 
energy distributions (SED) with a physically-motivated dust model to access the dust properties: dust mass, total-IR luminosity, 
polycyclic aromatic hydrocarbon (PAH) mass fraction, dust temperature distribution, and dust-to-stellar mass ratio. 

Results. Using a different SED model (modified blackbody), different dust composition (amorphous carbon in lieu of graphite), or a 
different wavelength coverage at submm wavelengths results in differences in the dust mass estimate of a factor two to three, showing 
that this parameter is subject to non-negligible systematic modelling uncertainties. We find half as much dust with the amorphous 
carbon dust composition. For eight galaxies in our sample, we find a rather small excess at 500 pm (< 1.5 cr). We find that the dust SED 
of low-metallicity galaxies is broader and peaks at shorter wavelengths compared to more metal-rich systems, a sign of a clumpier 
medium in dwarf galaxies. The PAH mass fraction and the dust temperature distribution are found to be driven mostly by the specific 
star-formation rate, sS FR, with secondary effects from metallicity. The correlations between metallicity and dust mass or total-IR 
luminosity are direct consequences of the stellar mass-metallicity relation. The dust-to-stellar mass ratios of metal-rich sources follow 
the well-studied trend of decreasing ratio for decreasing sS FR. The relation is more complex for highly star-forming low-metallicity 
galaxies and depends on the chemical evolutionary stage of the source (i.e., gas-to-dust mass ratio). Dust growth processes in the ISM 
play a key role in the dust mass build-up with respect to the stellar content at high sS FR and low metallicity. 

Conclusions. We conclude that the evolution of the dust properties from metal-poor to metal-rich galaxies derives from a complex 
interplay between star-formation activity, stellar mass and metallicity. 

Key words. ISM:evolution- galaxies:dwarf - galaxies:evolution- infrared:ISM- ISM: dust,extinction 


1. Introduction 

The processes by which galaxies evolve from primordial envi¬ 
ronments to present-day galaxies are still widely debated, but the 
seeds of this evolution lie in the star-formation histories of the 
galaxies, and in their interaction with their environment through 


gas infall, outflows or mergers. The interstellar medium (ISM) 
plays a key role in this evolution, being the site of stellar birth 
and the repository of stellar ejecta. 

Although interstellar dust represents only ~ 1% of the to¬ 
tal mass of the ISM it is an important agent in star formation. 
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Dust absorbs the stellar radiation that would otherwise dissoci¬ 
ate molecules and thus participates actively in the cooling of the 
ISM. It is also a catalyser for molecular gas formation by provid¬ 
ing a surface where atoms can react (Hasegawa & Herbst| 


[Vidali et al.|2004HLe Bourlot et al.|2()12[|Bron et al.|2014[ ). The 

presence of dust can increase the H 2 formation rate by about two 
orders of magnitudes compared to H 2 formation without dust 
(Tielens|2005 1 , thus facilitating star formation. 


Dust forms from the available heavy elements in the explo¬ 
sively ejected material from core-collapse supemovse (SN) and 


in the quiescent outhows from low-mass stars (Todini & Ferrara 

2001 

[Gomez et al.|2012a|bl 

Indebetouw et al. 120141 [Rowlands 

et al. 

2014[ Matsuura et al. 

2015J. The refractory dust grains 


coagulation in dense molecular clouds (Bazell & Dwek|ri990[ 


Stepnik et al. 2001 2003 [ Kohler et al.| 2012 2015| l, locking 


even more heavy elements in the solid phase of the ISM ( [Savage 
& Sembach|1996| Whittet|2003 1. Through destructive processes 


(such as erosion or sputtering, see Jones et al.|1994| 'l996||Sem 


[Dfaz-Cano & Jones|2008t [Bocchio et al.|2012||2014] l, elements 

are released again into the gas phase. Metallicity, defined as the 
mass fraction of heavy elements, or “metals”, in the ISM, is thus 
a key parameter in studying the evolution of galaxies. 

Understanding how dust properties evolve as a function of 
metal enrichment can provide important constraints for galaxy 
evolution studies. Dwarf galaxies in the local Universe are ideal 
targets for such a study as many of them have low metallicity and 
high star-formation activity. As such, they present star-formation 
properties and ISM conditions that are the closest analogues to 
those thought to be present in the primordial environments of the 
early Universe ( Madau et al. [19981 . 

Previous studies already demonstrated that the dust proper¬ 
ties in low-metallicity galaxies were notably different from that 
of metal-rich sources. Low-metallicity galaxies harbour warmer 
dust (e.g., Helou||1986 Melisse & Israel 1994[ [Galliano et’aT 


2005 [ Rosenberg et al. 2006[ Cannon et al.[[2066| l, lower poly 


cyclic aromatic hydrocarbon (PAH) abundances (e.g.. Madden 

2000 

Boselli et al. 2004| Engelbracht et al. 20051 |Wu et al. 

2006 

[O’Halloran et al.|20061 [Draine et al.|2007[ [Galliano et al. 

2008 

[Gordon et al. 2008 Wu et al.|201 l[[Sandstrom et al. 2012[l, 

higher gas-to-dust mass ratios (e.g., Issa et al.[[1990 Lisenfeld 

& Ferraral 19981 Hirashita et al.|2002l James et al.|20021 Draine 

et al.[2007[ Engelbracht et al.[2008 

Galliano et al.[2008[ [Munoz- 

Mateos et al.| 

20091 [Bendo et al.[ 

20101 [Galametz et al.|[2011[ 

Magrini et al. 

20111). Submillimetre (submm) excess emission. 


presently not accounted for by dust models, is observed in nu- 


2003112005 

Dumke et al.||20041 

Bendo et al.| 20061 Galametz 

et al. 12009 

2012 

|Zhu et al. 12000 

Bot et al. 120101 Gross! et al. 


Galametz et al.[[2()T4 Gordon et al. [[2014^1. Low-mass galaxies 


also show a broadening of the infrared (IR) peak of the SED 
(e.g., [Boselli et al.|2012[ Smith et al.|2012a[ Ciesla et al.|2014 , 


and a flattening of the far-infrared (FIR) slope (e.g., Boselli et a 
[2010[[^rtese et al.|2014[ l. However, all of these studies do not 
extend to very low-metallicities and/or have a limited number of 
sources below ~l/5 Zgp](12-i-log(0/H) = 8.0). 

In [Remy-Ruyer et al. ( 2013[ ), we derived the dust proper¬ 
ties in a systematic way for a large number of galaxies (109), 
among which more than half are dwarf galaxies, and ~ 35 % 
have Z < 1/5 Zq. We confirmed on this significant sample of 


* Throughout the paper, we assume 
12-Hlog(0/H)o= 8.69 [Asplund et al. 2009 1 


(0/H)e 


- 4.90xl0\ 


i.e.. 


low-metallicity sources that the dust is warmer at low metal- 
licities and identified several galaxies with submm excess. This 
study was made using only Herschel data and a modified black- 
body to model the dust emission. In [Remy-Ruyer et'd? (2014a 1 , 
we confirmed the higher gas-to-dust mass ratios (G/D) at low 
metallicities using the sample of Remy-Ruyer et al. ( 2013| l, and 
a semi-empirical spectral energy distribution (SED) model over 
the whole IR range. We showed that the G/D is actually higher at 
low metallicities than that expected from a simple description of 
the dust evolution in the ISM. The large scatter in the observed 
G/D is intrinsic to the galaxies and reflects a non-universal dust- 
to-metal mass ratio (s ee also[Dwek 1998[ De Cia et aT[[2013[ 


[Zafar & Watson|2013[ [Mattsson et al.[[20i4a| l. We showed that 

the metallicity was the main driver of the G/D but that the scat¬ 
ter was controlled by the different star-formation histories of the 
sources. Thus metallicity is not the only parameter shaping the 
dust properties. 


As a follow-up of Remy-Ruyer et'aL| ( [2013] 2014a 1 , we aim 
at providing a consistent picture of the evolution of the dust 
properties from metal-poor to metal-rich galaxies. On the obser¬ 
vational side, we present the complete catalog of IR-to-submm 
flux densities for the Dwarf Galaxy Survey (DGS). We use a 
semi-empirical dust SED model that accounts for starlight in¬ 
tensity mixing in the ISM, to interpret the whole IR-to-submm 
observed SEDs. In this work, we extend the range of dust prop¬ 
erties and look at the dust mass, the total infrared (TIR) luminos¬ 
ity, the PAH mass fraction and the dust temperature distribution. 
The dust temperature distribution is directly linked to the SED 
shape and will provide valuable insight on the average dust tem¬ 
perature, but also on the broadness of the FIR peak of the SED. 
We also discuss the evolution of dust-to-stellar mass ratios and 


link it to the evolution of the G/D observed in Remy-Ruyer et al. 
( 2014a| l. For the analysis we do not consider only metallicity as 
our main parameter, but also include the specihc star-formation 
rate and the stellar mass. 

The paper is organised as follows: in Section]^ we present 
the sample used for the study and the whole IR-to-submm 
dataset. Sectionj^presents the dust model, together with the free 
parameters (Sects. [TT]and[3:^, and errors on the best-ht param¬ 
eters (Sect. [33] l. In this section we also address how the use of 
a realistic dust model (Sect. [3.4[ |, or a different interstellar radia¬ 
tion held (ISRF, Sect.[3.5[), or a different dust composition (Sect. 


3 . 61 , or a different wavelength coverage at 500 //m (Sect. 3.7 1 


impacts the estimated dust properties (especially the dust mass). 
Section]^ is dedicated to the derived dust properties and their 
variation with metallicity and star-formation activity. In Section 
1^ we discuss our results on the dust-to-stellar mass ratio in the 
context of galaxy evolution. 

The paper also contains several appendice^in order to facil¬ 
itate the reading of the main body of the paper. Interested read¬ 
ers can hnd more details there. Appendix[A|presents how we up¬ 
dated Herschel photometry since the studies of [Dale et al.[ ( [2012| 
and [Remy-Ruyer et al.[ ( [2013[ l. Observing logs are given in 
Appendix [B[ for the DGS IRAC and IRS Spitzer data. Appendix 
[^compares IRAC and WISE photometry for the DGS sources to 
other measurements in the literature. Appendix[D|explains in de¬ 
tail the data reduction and preparation of the IRS spectra. Special 
cases of SED modelling are discussed in Appendix]^ 


^ online version only 
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A. Remy-Ruyer et al.: Linking dust emission to fundamental properties in galaxies 


2. Samples and Observations 

In this section we present the sample of galaxies and the set 
of IR-to-submm data used to build their observed dust SEDs. 
We also derive star-formation rates (SFR) and specific star- 
formation rates {sSFR) to quantify the star-formation activity. 


2.1. Samples 

For this study, we combine two samples of local galaxies ob¬ 
served with Herschel: the Dwarf Galaxy Survey (DGS, Madden| 
et al. 2013[) and the Key Insights on Nearby Galaxies: a Far- 


Infrared Survey with Herschel (KINGFISH, [Kennicutt et al.| 
201 l|l. The DGS is a sample of 48 star-forming low-metallicity 


local dwarf galaxies dedicated to the study of the dust and the gas 
properties in low-metallicity environments. The DGS sample in¬ 
cludes objects among the most metal-poor galaxies in the local 
Universe, with metallicities ranging from Z ~ 0.03 Zq to 0.55 Zq 


(Fig.[T]i. Stellar masses ra nge over 4 orders of magnitudes, from 
3x10° Mq to 3xl0'° Mq ( Madden et al.|2014 1 . 


The KINGFISH sample has been built from the Spitzer 
Infrared Nearby Galaxies Survey (SINGS, Kennicutt et al. 
|2003| l, and enables us to span a wider metallicity range by in¬ 
cluding more metal-rich galaxies. KINGFISH contains 61 galax¬ 
ies, including mostly spiral galaxies together with several early- 
type and dwarf galaxies. The metallicities in the KINGFISH 
sample range from Z ~ 0.07 Zq to 1.20 Zq. The two samples 
complement each other very well in terms of metallicities and 
stellar masses (Fig.[2l. 

The KINGFISH sources have been classified in terms of 
'‘star-forming” or “Active Galactic Nuclei” (AGN) type by 


Moustakas et al.|( 2010i using an optical emission-line diag 
nostic diagram (Baldwin et al. 19811. According to [Kennicutt 


et al. (201 l|l, the detected AGNs are low-luminosity AGNs 


and none of the galaxies have a dominant AGN, except for 
NGC 1316. In terms of mid-infrared (MIR) emission line di¬ 


agnostics, NGC 4725 has a [O iv]/[Ne ii] ratio of ~ 1 (Dale et al. 
|2006| l and can be considered to be an AGN-dominated source. 
NGC 4736 has [O iv]/[Ne ii] ~ 0.3 and is an “intermediate” ob¬ 
ject with a AGN contribution of ~ 50%. The other KINGFISH 
AGN sources have [O iv]/[Ne ii] <0.1 or no detection of the 
[O iv] 25.9 pm line, equivalent to an AGN fraction < 10 % 
(Genzel et al. 1998[ |Tommasin et al. |2010 i. Except for these 
three galaxies, we thus do not expect the AGN to significantly 
impact the IR emission in the KINGFISH galaxies. 

The metallicities in both samples have been determined in 
Madden et al. (2013|l and [Kennicutt et "aL (2011 1 using the 


strong-line” calibration from |Pilyugin & Thuan| ( |2005] l, linking 
a specific oxygen optical lines raticr to the oxygen abundance in 
terms of 12+log(0/H). We choose this specific metallicity cali¬ 
bration because it is available and homogeneous for both sam¬ 
ples, in order to avoid any biases by using different calibrations. 
The stellar masses were derived by Madden et al. ( |2014 1 for the 
DGS using the prescription of |Eskew et al.| ( |2012[ ) and the IRAC 
3.6 pm and 4.5 pm flux densities. We use this prescription for the 
KINGFISH sample for consistency purposes and to reduce the 
scatter in the dust-to-stellar mass ratios (as seen in |Remy-Ruyer 
jet al.|2014b| ). When no IRAC observed flux densities are avail¬ 
able we use synthetic flux densities computed from the best-fit 
SED model (see Section]^. Coordinates, distances, and metal¬ 
licities for the DGS can be found in|Madden et al.|(|20131l, and in 
Kennicutt et aT] ( |200^[20TT] ) for the KINGFISH sample. Stellar 


The R 23 ratio, R23=([OII]43727-t[OIII]444959,5007)/H/3. 


DGS, KINGFISH - Metallicities 



DGS, KINGFISH - Stellar masses 



Log[M,,,yMo] 


Fig. 1. Metallicity (top), and stellar masses (bottom) distribu¬ 
tions of the DGS (purple), KINGFISH (orange) samples. The to¬ 
tal distribution is indicated in grey. In the top panel, solar metal¬ 
licity is indicated here to guide the eye, as well as the 1/50, 1/10 
and 1/5 Zq values. 


masses are tabulated in Table We have a total of 109 galaxies 
for this study. 


2.2. Herschel data 


The 109 DGS and KINGFISH galaxies have been observed 
with Herschel in 6 photometric bands at 70, 100 and 160 pm 
with the Photodetector Array Camera and Spectrometer (PACS, 
Poglitsch et al.| 2010| l, and at 250, 350 and 500 pm with the 


Spectral and Photometric Imaging REceiver (SPIRE, Griffin 
let al.|2010D . 

The data reduction and photometry for the DGS and 
KINGFISH galaxies are presented in |Remy-Ruyer et'ST ( 2013| l 
and Dale et aL] ( |2012) l, respectively. However, the data reduction 
techniques and the calibration of the instruments have been up¬ 
dated since the publication of|Remy-Ruyer et ah (201 3|l and Dale 
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A. Remy-Ruyer et al.: Linking dust emission to fundamental properties in galaxies 


jet al.| (20 T^ . In this work we therefore present the latest updated 
data. For the DGS data, the Scanamorphos PACS maps are re¬ 
processed with version 23 of Scanamorphos ( Roussel|2013[ ), and 
the SPIRE maps are reprocessed with version 12 of the Herschel 
Processing Interactive Environment (HIPE - Ott|2010 l. We find 
an average difference of 5%, 8% and 2 % at 70, 100 and 160 
^m between the updated flux densities and the flux densities of 
Remy-Ruyer et al.](|2013|l, and in some extreme cases the differ¬ 


ence is as high as 30%. Eor the SPIRE flux densities, the average 
difference is 3%, 5% and 3% at 250, 350, and 500 yum respec¬ 
tively. The updated PACS and SPIRE flux densities for the DGS 
galaxies are presented in Table 

Eor the KINGFISH data, we update the flux densities pre¬ 
sented in Dale et al. ( 2012[ ), with an average difference with the 
|Dale et al.| ( |2012^ ux densities of 7% at 160 //m, and ~ 2% 
for SPIRE wavelengths. Detailed information regarding the up¬ 
date of the DGS and KINGEISH Herschel data can be found in 
Appendix Colour corrections are applied to the SPIRE flux 
densities in both samples to account for the dependence of the 
SPIRE beam areas on the FIR spectral shape of the sources (see 
Appendix [a). 


Handbook (version 2.0.1 Section 4.11.l|^ and depends on the 
source aperture radius. Following the recipe given in the IRAC 
Instrument Handbook, we do not apply this correction to small 
and compact sources (see Table|^. 


The uncertainty on the flux density is computed by summing 
in quadrature the calibration error, the error from the background 
determination and the error from the source flux determination. 
We adopt a calibration error for the four IRAC bands of 10% as 
recommended in the IRAC Instrument Handbooljf] We consider 
that galaxies are not detected when the computed flux density is 
lower than three times its corresponding uncertainty in a given 
band. The final 3cr upper limit is reported in Table 


IRAC photometry is available in the literature for 29 
DGS sources; Hunt et al. ( 2006|l; Dale et akl ( |2007[ |2009a| l; 
Engelbracht et al. ( 2008) 1; Galametz et al. ( 2009| . If we compare 


our measurements to those from the literature we get a fairly 
good agreement between them, with some outliers. The details 
of the comparison with the literature measurements and possible 
explanations for the outliers are given in Appendix [C| 


2.3.2. IRS 


2.3. Spitzer data 


We collect Spitzer observations for our sample in order to con¬ 
strain the warm dust and the PAH emission in MIR wavelengths, 
with photometric data from the InfraRed Array Camera (IRAC), 
and from the Multiband Imaging Photometer for Spitzer (MIPS), 
and spectroscopic data from the InfraRed Spectrograph (IRS). 
Spitzer observations of several DGS sources were already avail¬ 
able in the Spitzer databas^ New complementary Spitzer obser¬ 
vations for 19 DGS galaxies have also been obtained during the 
cycle 5 program to complete the set of existing Spitzer data (Dust 
Evolution in Low-Metallicity Environments: PI. F. Galliano; ID: 
50550), and for which only the MIPS data at 24, 70 and 160 
yum have been published so far (in |Bendo et al.|2012| l. Thus we 
only present the IRAC and IRS observations for the DGS sample 
(Observing Logs in Appendix [B]l. 

Spitzer IRAC and MIPS photometry for the KINGFISH 
galaxies is taken from |Dale et al.| ( 2007| l. We do not attempt 
to use the IRS data existing for the KINGFISH galaxies, as 
the galaxies are very extended and a complete coverage of the 
sources with IRS is not available (|Dale et al.|2009b|l. 


2.3.1. IRAC 

IRAC data at 3.4, 4.5, 5.8 and 8.0 yum is available for the 
DGS galaxies, and the IRAC maps are retrieved from the 
NASA/IPAC ISA database for Spitzer data. For three galax¬ 
ies, HS 1442H-4250, SBS 1415-H437 and UM 311, only two 
wavelengths are available (4.5 and 8.0 pm). After subtracting 
the background and possible contamination from background 
sources or foreground stars, flux densities are extracted from the 
maps using aperture photometry. In most cases, we use the same 
apertures as those used for Herschel photometry. The final IRAC 
flux densities and the apertures used are given in Table 

The IRAC calibration is based on point-source photome¬ 
try for a 12" radius aperture. An additional aperture correc¬ 
tion is needed to account for the emission from the wings of 
the PSF and the scattering of diffuse emission across the IRAC 
focal plane. This correction is given in the IRAC Instrument 


The query form is available at: http://sha.ipac.caltech.edu 


The DGS galaxies have a very rapidly rising continuum in the 
MIR, and the IRS spectra gives important constraints on the con¬ 
tinuum from 20 to 40 pm. Moreover, the spectra can also be used 
to put constraints on the PAH emission (see Section 
intensities have been extracted by Cormier et al.| (j^ 


3.2 1 . Line 
5]ifor the 


[SIV] 10.5 pm, [Nell] 12.8 pm, [Nelll] 15.6 yum, [SIII] 18.7 pm 
and 33.5 yum spectral lines, in the compact sources of the DGS. 

The IRS spectra for the DGS galaxies have been extracted 
from the Cornell AtlaS of Spitzer Infrared spectrograph Source^ 
(CASSIS v5, Lebouteiller et al.|201 l| l. Most of the DGS galaxies 
were observed using the staring mode, except for a few extended 
sources observed with the mapping mode, for which the spec¬ 
tra were reduced manually with CUBISM ( Smith et al.|2007a i. 
We were able to obtain Short-Low (i.e., short wavelength data 
at low spectral resolution, SL) and Long-Low (i.e., long wave¬ 
length data at low spectral resolution, LL) data for 43 galaxies 
in total. More details about the extraction and the data reduction 
of the IRS spectra can be found in Appendix [P] 

It is then necessary to rescale the SL and LL spectra in order 
to match the photometry. We derive synthetic IRS photometry to 
correct the two modules and use all of the constraints we have 
in this wavelength range; the IRAC 5.8 and 8.0 yum, WISE 12 
and 22 pm (see Section [274) 1 and MIPS 24 pm bands. IRAC 5.8, 
8.0 and WISE 12 pm are used simultaneously to derive a correc¬ 
tion factor for the SL module that depends on the wavelength. 
However, for the LL module, the two constraints do not sample 
well the LL spectrum, and the LL correction factor is thus a con¬ 
stant. We use MIPS 24 pm in most cases (or WISE 22 pm when 
not available). SL and LL can be treated separately as they are 
two independent observations. More details about the rescaling 
of the IRS spectrum can be found in Appendix [^ The median 
correction applied at 5.8 pm, 8.0 yum, 12 yum, and 24 yum are 
1.7, 2.0, 1.0 and 1.2 respectively. Note that this correction as¬ 
sumes that the spectral shape of the area observed by the IRS 
slits likewise describes the expected spectral shape of the full 
galaxy. This is true for compact sources but can be erroneous for 


^ This document is available at http://irsa.ipac.caltech.edu/data/ 
SPITZER/docs/irac/iracinstrumenthandbook. 

® This can be found in Section 4.3 of the IRAC Instrument Handbook. 
^ The database is available at http://cassis.astro.comell.edu/atlas/ 
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more extended galaxies, except if the region falling within the 
IRS slits dominates the total emission of the galaxy in the MIR. 

To improve the quality of the noisiest spectra (e.g. 
HS 1304+3529, Fig |D.l| top panel), we smooth the spectra until 
we reach a signal-to-noise ratio (S/N) of 3 at every wavelength. 
This smoothing step is applied for 22 DGS galaxies. However, 
for three galaxies, HS 1236+3937, HS 1442+4250, Tol 0618- 
402, even with the smoothing step, we cannot reach a S/N of 
3 for any point in the spectrum. Thus we do not consider these 
spectra in the SED modelling. In two galaxies, HS 2352+2733 
and UGCA 20, the IRS slits are not centred on the source posi¬ 
tion, thus we do not present these spectra either. The remaining 
38 IRS spectra are shown in Appendix [P] Fig.|D.3| 


2.4. WISE data 


The Wide-held Infrared Survey Explorer (WISE, |Wright et al. 
|2010| l observed the whole sky at 3.4, 4.6, 12 and 22 jum, and 
gives an additional constraint at 12 fim, very valuable for probing 
the MIR range of the observed SED. The InfraRed Astronomical 
Satellite (IRAS, Neugebauer et al.||1984| l also provides a con¬ 
straint at 12 fim (see Section |2.5| l but was not able to detect the 
faintest dwarf galaxies. We also need the 12 and 22 /im con¬ 
straints to match the IRS spectra to the photometry (Section 
2.3.2|l. WISE Hux densities at 3.4 and 4.6 /rm are also given for 


completeness. We present here WISE data for the DGS galaxies 
only as Spitzer and IRAS already provide the equivalent spectral 
coverage for the brighter KINGEISH galaxies. 

The WISE maps were retrieved from the NASA/IPAC ISA 
database and the AllWISE databas43 Most sources in the DGS 
are resolved by WISE as the resolution is ~ 6" for the hrst three 
bands and 12" for the WISE 22 pm band. Prohle-ht and stan¬ 
dard aperture photometry provided by the WISE database can 
underestimate the brightness of resolved sources and can suffer 
from confusion with nearby objects (Section 2.2 of the AllWISE 
explanatory supplemenj^. 

We therefore perform our own aperture photometry on the 
maps following the method outlined in Section 4.3 of the All 
WISE explanatory supplement. The apertures are the same as 
for IRAC and Herschel photometry whenever possible. In some 
cases, the morphology of the source in the hrst two bands is 
quite different from the morphology of the source in the last two 
bands. These two sets of wavelengths are not tracing the same 
physical component, old stellar population at 3.4 pm and 4.6 pm 
on one side, warm dust at 12 pm and 22 pm on the other side, and 
the near-infrared (NIR) morphology may differ from the MIR 
morphology. In this case we have to use two different apertures 
if we want to encompass the whole emission of the source and 
optimise the S/N ratio at all four wavelengths. Eor the unresolved 
sources, we use the prohle-ht photometry provided by the All 
WISE database. 

The uncertainty comes from the hux determination, from the 
background estimation, from correlated noise in the maps, and 
from the calibration and measurement of the zero-point magni¬ 
tudes (see the AllWISE explanatory supplement). The calibra¬ 
tion of WISE have been tied to that of Spitzer IRAC 3.6 pm for 
WISEl, IRAC 4.5 jum for WISE2, IRS f or WISE3, and MIPS 
24 pm for WISE4 (see Jarrett et al.|201 l| l, thus we adopt a cal¬ 
ibration error of 10 % at 3.4 pm and 4.6 pm, 5% at 12 pm, and 
7% at 22 pm. As for IRAC photometry, we report upper limits 


* The query form is available at: http://sha.ipac.caltech.edu 

* This document is available at http://wise2.ipac.caltech.edu/docs/ 
release/allwise/expsup/index.html. 


if the detection is below 3cr. The hnal WISE hux densities and 
apertures used are given in Table 

WISE measurements compare well with other instruments, 
IRAC at 3.6 pm and 4.5 pm, IRAS at 12 pm, and MIPS at 24 
pm. The details of the comparison are given in Appendix [C] 


2.5. All-sky survey data from 2MASS and IRAS 

We complete our set of Spitzer WISE and Herschel data by 
searching the literature for data from 2MASS (Two Microns All 
Sky Survey, Skrutskie et al.|2006| and IRAS. 

The 2MASS data in the J (1.24 /tm), H (1.66 pm) and Kj 
(2.16 pro) bands, for the DGS sample is given in Table and 
has been compiled from the literature: the NASA/IPAC ISA 
2MASS Point Source CataloJ^ the 2MASS Extended Objects 
Einal Release, the 2MASS Large Galaxy Atlas (Jarrett et al. 
20031, Engelbracht et al. ( |2008| l; [Dale et ah P009a I. In some 
cases, the original data is given in magnitudes. To convert this 
into Hux densities, we use the zero-magnitude Hux values from 
Cohen et al.|(p00^. 2MASS d ata for the KINGFISH galaxies is 
presented in [Dale et al.| ( |2007| . 

We also compile IRAS data for the DGS at 12, 25, 60 and 
100 pm from the literature; the NASA/IPAC ISA IRAS Faint 
Source (v2.0 1990)® and Point Source (v2.1) Catalogs® , |Rice| 
et al. ( |1988| l; [Sanders et al. (2003 i; Engelbracht et al. ( [2008 1 
and is given in Table [^ Note that IRAS is not very sensitive for 
the low-luminosity dwarf galaxies and that the resolution varies 
from ~ 30" at 12 pm to ~ 120"at 100 pm. Even for the detected 
sources, such low resolution implies that several sources may be 
mixed in the beam and indistinguishable. That is why some of 
the IRAS huxes of Table have been noted as unreliable (with 
the note ^). 25 DGS galaxies have IRAS measurements and only 
15 are detected in all of the IRAS bands. KINGEISH IRAS data 
have been extracted from [Sanders et al.j P003] l and the IRAS 
Eaint Source (v2.0 1990) and Point Source (v2.1) Catalogs. 


2.6. Star formation activity 
2.6.1. Star-formation rates 


To estimate the SER for our sources, we consider two of the 
most widely used SER tracers, the far-ultraviolet (EUV) and the 
Ha luminosities. The observed stellar emission traced either by 
EUV or Ha has to be corrected for dust attenuation, and this 
correction is usually made using observational dust tracers such 
as the TIR or 24 pm lumino sities ([Kennicutt et al.[2009[|Calzetti| 
[et al.|20f0l|Hao et al.|201l1|Kennicutt & Evans|2()12[). 

In our sample, we can derive SER for 90% of the galaxies us¬ 
ing Ha+TIR, 88% using Ha+24 pm, 77% using EUV+TIR and 
78% using EUV+24 pm. We thus use the diagnostic based on the 
observed Ha luminosity, corrected for attenuation using Cjir. 
and the formula by [Kennicutt et al.[ ([2009|l for a Kroupa initial 
mass function (IMF). Moreover, [Kennicutt et al.[ ( |2b09| l showed 
for the Spitzer-SINGS sample that this composite tracer gave the 
most robust SFR measurements compared to Ha measurements 
corrected for dust attenuation using the B aimer decrement. 

The integrated Ha luminosities for the sample have been 


taken from the literature (mostly from Moustakas & Kennicutt 
|2006[ [Kennicutt et al.||2008[ |2009| ), and are listed in Table [4( 
Care has been taken in providing Ha luminosities corrected for 
underlying stellar absorption, for Nil line contamination and 
foreground Galactic extinction. The same extinction curve was 


The catalog is available at: http://irsa.ipac.caltech.edu 
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used to correct our sources for foreground Galactic extinction: 
A(Ha)=0.6 Ab ( |0’DonneII||1994[ ). The B-band extinction. An, 
is taken from the NASA/IPAC Extragalactic Database (NEEp^ 
and the measurements of |Schlafly & Einkbeiner| ( |2011| l. Eor the 
hve galaxies missing either Ha or TIR data, we use the mean of 
the other SER estimates. The error on the SER is derived by tak¬ 
ing the standard deviation of all of the SER estimates when more 
than one is available, otherwise we use the 20% median error on 
the SER in the sample (reported in Table [T]). 

|Lee et all] ( |2009| l cautioned the use of Ha to derive SER for 
galaxies with integrated Lna < 2.5 X 10^® erg.s ' or equivalently 
SER < 0.01 Mq yr '. Ten galaxies in our sample are below these 
limits in Lna or SER. Eor these very faint galaxies, Ha-SER 
will be underestimated, even after correction for dust attenua¬ 
tion ( |Lee et ah 2009 Meurer et al.|2009|. Eor th ese ten galaxies, 
we use the calibration provided by Lee et al. (20091 (specifi¬ 
cally their equation 10), which converts the non-dust corrected 
Ha luminosities into dust-corrected SER, based on the |Kennicutt| 
( |1998| l EUV-SER. These obtained SERs are from a calibration 
assuming a Salpeter IME ( Kennicutt|1998 1 and we thus convert 
them into Kroupa IME-based SERs using the correction factors 
tabulated in Table 1 of|Kennicutt & Evans |(|2012|). The error on 


the SER for these ten galaxies is taken from Table 2 in Lee et al. 
(2009|l. The resulting SERs for these ten galaxies are a factor of 


2 larger than the previously estimated SERs from Ha and TIR. 
The Ha luminosities and hnal SER are given in Table and 
shown in Eig.|^ 

The SERs in our sample cover four orders of magnitude, with 
a median SER of 0.27 Mgyr ' and 0.54 Moyr ' for the DGS 
and KINGEISH samples respectively (see Table[^. We see from 
Eig.j^that at a given stellar mass, galaxies with the highest SERs 
have the lowest metallicities. This is consistent with the hndings 
of|Lara-L6pez et aT ( 2010| l; Mannucci et al. ( 2010[ ); [Hunt et al. 


(2()12|l on the “fundamental metallicity relation” (EMR) describ 


ing the various scaling relations between metallicity, stellar mass 
and star-formation rate. 


2.6.2. Specific star-formation rate 

We also consider the specific star-formation rate, sSFR, i.e., 
SER normalised by the stellar mass, to remove scaling effects in 
our analysis of such a wide variety of environments. The sSFR 
values are given in Table and plotted as a function of metallic¬ 
ity in Eig.|^ 

The median sS FR is one order of magnitude higher in the 
DGS than in KINGEISH (0.51 Gyr"' versus 0.04 Gyr '). In 
the KINGEISH sample, we distinguish the galaxies according 
to their nuclear types as dehned by |Kennicutt et al. (201 l)l: 


star-forming (SE) versus non-thermal emission (AGN). Th 
KINGEISH “AGN” galaxies have the lowest sSFR in Eig. ^ 
with a median sSFR of 0.01 Gyr \ versus 0.06 Gyr“^ for the 
KINGEISH “SE” sample. The stellar masses were estimated on 
the IRAC 3.6 fj.m and 4.5 /urn flux densities, which might not 
be appropriate for AGN-type galaxies, possibly overestimating 
Mstai and thus yielding underestimated sSFR values. 

The clear correlation between sSFR and metallicity seen 
in Eig. 1^ follows directly from the EMR. The sS FR can also 
be interpreted as the amount of star formation occurring today 
compared to the output of the cumulated past star-formation, 
i.e, the mass of stars already formed, tracing the integrated star- 
formation history. The sS FR value can thus be boosted up by 
the low stellar mass of the dwarf galaxies even if the present 


star-formation rate is comparable to that in more metal-rich en¬ 
vironments. This is consistent with the so-called “downsizing” 
effect, i.e., more massive galaxies are more efficient in convert- 


Cowie et al. 1996 [Gavazzi & Scodeggio 

1996| Boselli et al. 

2001 

Gavazzi et al.|2002| Brinchmann et a 

.|2004| Hughes et al. 


20131. The less massive galaxies form the bulk of their stel¬ 
lar mass later than more massive systems, and thus undergo a 
slower metal enrichment. This is consistent with our picture of 
low-mass, low-metallicity galaxies having a higher sSFR today 
than their massive metal-rich counterparts. 


3. Dust modelling in the DGS and KINGFISH 
galaxies 

In this section, we describe the dust SED model, the parameters 
and the derivation of the errors on the best-ht parameters (Sects. 


3.2 3.3 I. In Sect.|3.4| we compare our results with the dust 


properties derived with a modified blackbody model (namely the 
dust mass and the temperature). In Sect. 3.5 and |3.6| we inves¬ 
tigate the effect of a different description for the radiation field 
and of a different dust composition on the dust properties, re¬ 
spectively. And hnally in Sect. 3.7 we show that the 500 /rm data 


point is a necessary constraint to accurately estimate the dust 
mass, and investigate the influence of a potential submm excess. 

3.1. Description of the SED modei 

We adopt the semi-empirical dust SED model presented in 
Galliano et al. ( |201 l| l. The model adopts the Galactic grain 
composition made of silicate grains, carbon grains in the form 
of graphite, and PAHs. The optical properties are taken from 
Weingartner & Draine ( |2001| l; Laor & Draine (1993 i and Draine 
& Li (2007|l respectively. The assumed size distribution is that 


determined by |Zubko et al.| ( [2004 1 for their BARE-GR-S model, 
but we allow the relative normalisation of the PAH component 
to vary for more flexibility. The graphite-to-silicate ratio is kept 
fixed to the [Zubko et al.| ( |2004[ ) value. We assume that the ISRE 
illuminating the dust grains has the spectral shape of the so¬ 
lar neighbourhood ISRE ( Mathis et al.ll983 i. Only its intensity 
varies, controlled by the parameter The temperature fluc¬ 
tuations of stochastically heated grains are computed using the 
transition matrix method ( [Guhathakurta & Draine|1989] ). 

We assume that the dust properties are uniform, i.e., the size 
distribution and the grain optical properties are constant within 
the galaxy, and that only the starlight intensity varies. We dis¬ 
tribute the mass into different mass-elements of uniform illu¬ 


mination with the empirical recommendation from Dale et al. 
(2001 1 : the distribution of dust mass per starlight intensity can 
be approximated by a power law of index a, dM/dU oc 
with U varying between U^nin and t/min+ At/. This simply ex¬ 
presses that most of the dust mass should reside in the coldest 
components, with low starlight intensities U. In most cases, this 
is flexible enough to reproduce dense and diffuse media, and pro¬ 
vides a simple parameterization of the physical conditions in the 
ISM. 

Emission from old stars can also contribute to the IR emis¬ 
sion and especially in the NIR. Thus we add a stellar continuum 
to the dust emission, parameterised by the stellar mass of the 
galaxy M*. This parameter is not designed to estimate the stellar 


This database is available at http://ned.ipac.caltech.edu/ 


U is expressed relative to the solar neighbourhood value Uq = 2.2 
X 10^^ W.m-2 
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Fig. 2. Star-formation rates, SFR, (left), specific star-formation rates, sS FR, (right), for the DGS (crosses) and KINGFISH samples 
as a function of metallicity, colour coded with Mjtar (colour legend in right panel). For the KINGFISH sample, we distinguish 
the sources that have been classified as “star-forming” (SF, downward triangles) or as “AGN” (stars) by |Kennicutt et aL ( 201 l| l. 
The distribution of each parameter is indicated on the side of each panel for both samples: solid line for DGS and dashed line for 
KINGFISH. 


mass, but rather to have a good fit in the NIR (see Galliano et al. 
|2011| for details). 

For each band at a wavelength T, the synthetic luminosity is 
computed by convolving the model with the spectral response of 
each band and using the appropriate spectral convention of each 
instrument. This is done for all of the observational constraints 
except at SPIRE wavelengths, as the input flux densities are al¬ 
ready colour-corrected (see Sect. |2.2[ ). 

The fit is performed with the mpcurvefit iterative procedure 
in IDL based on the Levenberg-Marquardt method. This model 
has been used to model galaxies before by|Galametz et al.|P0091 
20T0l[20T7]|2013a|l;|Hony et al.| ( p0T0| l; [Meixner et al.| ( |20101 l; 


O’Halloran et al.| ( |2010| l. 


The Galactic PAH mass fraction is /pah0 =4.57% from 


Zubko et al. 


IRS spectra to constrain it. These two parameters are used in 
order to allow for a more flexible and a better fit in the MIR, es¬ 
pecially to reproduce the steeply rising continuum of the dwarf 
galaxies, but will not be discussed any further. 

From the final best-fit models, we derive the TIR luminos¬ 
ity, Ltir, by integrating the total dust SED between 1 and 1000 
fim for each galaxy. We compute as well the first and second 
moments of the starlight intensity distribution, (U) and (crlTf, 
corresponding to the mass-averaged starlight intensity and the 
variance in the starlight intensity distribution re spectively. (U) 
and (crUf are given by (see Eqs. 9 to 12 from 

[2OTT] ): 


Galliano et al. 


3.2. Setting the parameters 

The model can be described by six main free parameters: 

- Mdnst- the total dust mass, 

- Umin- the minimum of the starlight intensity distribution, 

- At/: the difference between the maximum and minimum of 
the starlight intensity distribution, 

- a: the index of the power law describing the starlight inten¬ 
sity distribution, 

- /pah: the PAH mass fraction, normalised to the Galactic PAH 
mass fractiorp^ and 

- M*: the stellar mass. 

Units and bounds used for the fitting are described in 
Galliano et al. ( 201 l| l. In addition, we vary the mass fraction 
of very small grains (non-PAH grains with radius a < 10 nm), 
/vsg, relative to the Galactic valu^^ It is constrained by the 8 - 
30 jum photometry, plus IRS spectroscopy for the DGS sources, 
and controls the MIR continuum. We also vary the ionised-to- 
total PAH mass ratio in the DGS sources, as we have the 


1 dMdust 

(U) =- Ux^^dU 

ALdu.st J Umin 


(cruy 


-^(dust Ju 


t/niin+AC/ 


dU 

{U-{U})^x^dU 
dU 


( 1 ) 


( 2 ) 


( |2004) . 

We compute the Galactic mass fraction of non-PAH grains with 
radius <10 nm in Zubko et al. 1 2004 1 model: /vsgo = 16%. 


(U) can be seen as an indirect measure of the average dust 
temperature and crU of the broadness of the FIR peak of the 
SED. 

This model is applied to the observed SEDs of DGS and 
KINGEISH galaxies. An example of fit is given in Eig.|^ the oth¬ 
ers are shown in Eigs.[^and[T^ For four DGS galaxies, there 
are not enough constraints to properly fit a SED: HS 1236-1-3937, 
HS 2352-P2733, Tol 0618-402 and UGCA 20. We do not con¬ 
sider KINGFISH galaxies without Herschel detections either: 
NGC 0584, NGC 1404, DDO 154 and DDO 165. We also re¬ 
move from the subsequent analysis the three KINGFISH galax¬ 
ies for which the AGN emission is dominant: NGC 1316, 
NGC 4725 and NGC 4736 (see Section[2T]i. Indeed, [Ciesla et ah] 
(2015 1 showed that the emission of an AGN can significantly im¬ 
pact the total IR emission from contribution >40%. For the other 
KINGFISH “AGN” sources, we do not expect the total IR emis¬ 
sion to be affected by the AGN contribution on global galaxy 
scales, and due to the difficulty of accurately constraining such 
low AGN fractions ( [Ciesla et al.|2015 1 we do not apply any spe¬ 
cific modelling. We have a total of 98 galaxies to which our SED 
model is applied. 
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NGC4449 



1 10 100 1000 
Wavelength [/^m] 


Fig. 3. SEDs for a DGS source NGC 4449. The observed SED 
includes the Herschel data (purple crosses) as well as any avail¬ 
able ancillary data (in orange). The different symbols code for 
the different instruments; Xs for 2MASS bands, stars for Spitzer 
IRAC and MIPS, diamonds for WISE and triangles for IRAS. 
The IRS spectrum is displayed in orange. The total modelled 
SED in black is the sum of the stellar (green) and dust (red) 
contributions. The modelled points in the different bands are the 
hlled blue circles. Residuals are shown on the bottom panel. 


3.2.1. Wavelength coverage 


The wavelength coverage is not entirely the same for all of the 
galaxies. Eor the DGS galaxies, the observational constraints 
come from the 2MASS survey and the Spitzer, WISE, IRAS 
and Herschel instruments (Section]^. Eor the KINGEISH galax¬ 
ies, the observational constraints originate from the same in¬ 
struments (except WISE and Spffzer/IRS). We have complete 
spectral coverage for 93% of the sources in the 1-5 pm NIR 
range (five wavelength^^, 87% in the 5-50 pm MIR range 
(four wavelengths), 90% in the 50 - 200 pm EIR range (four 
wavelengths), and 76% in the 200 - 500 pm submm range (three 
wavelengths). [Ciesla et al. ( 2014| l showed that SPIRE constraints 
were particularly important to account for the cold dust (see also 
Section [T7] l. The 24% of our sample without SPIRE detections 
are all dwarf galaxies, harbouring particularly warm dust, with 
a peak of the EIR SED at wavelengths ~ 30 - 70 pm, well sam¬ 
pled by constraints until 160 pm. Thus, although heterogeneous, 
the wavelength coverage achieved in our sample is excellent and 
we are conhdent in the parameters we derive with these sets of 
constraints. 


When several observations are available at the same wave¬ 
length we favour the fluxes obtained with higher resolution data. 
When IRS constraints are available, we weight the IRS data 
points for them to contribute equally as the other MIR con¬ 
straints in the fit. Whenever PAH features are absent from the 
IRS spectrum, we fix /pah = 0, and fix /,,„ = 0.5, to reduce the 
number of free parameters. Some galaxies are not detected at 
one or several wavelengths, and we impose to the best-fit model 
to be consistent with the upper limits. 


We consider only one datapoint when several observations are 
available for one wavelength (e.g., either WISE 3.4 pm or IRAC 3.6 
pm, either IRAS 100 pm or PACS 100 pm, etc.) to avoid redundancy. 


3.2.2. Additional features 


For eleven DGS galaxies the observed MIR continuum shape 
cannot be well reproduced by our model. An example is given 
in the top panel of Fig. for SBS 1533H-574. For these eleven 
sources, we add to the model an extra modified blackbody 
(modBB) component at MIR wavelengths, with a fixed /? = 2.0 
and a temperature varying between 80 and 300 K. It can be 
physically interpreted as a contribution from hot Hit regions to 
the total emission of the galaxy. As dwarf galaxies have small 
physical sizes and low dust attenuation, the emission from the 
energetic Hit regions can indeed affect the total emission from 
the whole galaxy (|da Cunha et al.||2008 Galametz et al.|[20T0 


[Hermelo et al.|2013| l. This is observed also for FIR-line emission 
that is often dominated on galaxy-wide scales by the emission 
from star-forming regions or diffuse ionised gas in dwarf galax- 
ies (e. g., Cormier et al. 2015| l. In these cases, the |Dale et ah 
(2001 1 empirical recipe of dMIdU oc f/~“ may not be flexible 
enough to provide a satisfactory fit. Of course, a modBB is not 
the optimal way of modelling the emission from hot grains in Hit 
regions (see |Galliano et al.|[2008] [Groves et al.||200^ as these 
grains are likely not in thermal equilibrium, but we adopt this 
approach as a first approximation. 



1 10 100 1000 


Wavelength [/^m] 



1 10 100 1000 
Wavelength [/.tm] 


Fig. 4. SED for SBS1533H-574 without (top) and with (bottom) 
the extra MIR modified blackbody, shown in yellow. The other 
colours and the symbols are the same as for Figj^ 
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This additional component does not impact significantly the 
total dust mass determination: for each of the eleven cases, dust 
mass estimates with and without the extra MIR modBB are con¬ 
sistent within errors except for SBS 1533-1-574. However, in this 
case the addition of a warm modBB dramatically changes the 
modelled MIR-to-submm shape of the SED (see Fig. 1^ The 
best fit MIR modBB temperatures are reported in Table]^ The 
fits without the MIR modBB are shown in Appendix]^ Fig. |E.l| 
for comparison. 

Specific modelling was required for five DGS galaxies and 
is detailed in Appendix]^ 


3.3. Estimating the errors on the best-fit parameters 

The errors on the best-fit parameters are estimated by generating 
300 random realisations of the observed SEDs within the obser¬ 
vational errors, taking special care for errors which are corre¬ 
lated between different bands (see below). For each galaxy, we 
perform fits of the 300 random realisations of the SEDs in order 
to obtain a distribution for the parameters. Following the method 
outlined in Remy-Ruyer et al. (20131, the perturbation to the ob¬ 
served fluxes is the sum of two components: 


- A normal random independent variable representing the 
measurement errors at each wavelength. 

- A normal random variable describing the calibration errors 
that takes into account the correlation between the wave¬ 
bands for each instrument (as detailed below). 


The errors are taken as the 66.67% confidence level. The 
best-fit parameters with errors are given in Table for all of the 
galaxies. 

We detail below the decomposition of the calibration errors 
and eventual correlations between the different bands for all 
of the instruments, except Herschel as this decomposition has 
already been presented in Remy-Ruyer et al. ( |2013[ ). Taking into 
account the correlations in the observationnal errors allows for 
smaller error bars on the dust parameters. 


2MASS: Jarrett et al. (2003 1 quote a 2 - 3% uncertainty on the 
zero-magnitude flux values. To be conservative, we assume an 
independent error of 3% in each band. 


I RAC: The total calibration uncertainty used for the DGS and 
KINGFISHIRAC fluxes is ~ 10%. This can be decomposed into 
two parts: 


Reach et al. (2005|l give a 2% uncertainty in all of the IRAC 


bands. This error is correlated between the four bands. 

The IRAC Instrument Handbook (Section 4.3) also gives a 
10% error to account for several systematic effects in the cal¬ 
ibration. This error is independent from band to band. 


IRS: The IRS spectra extracted from the CASSIS database pro¬ 
vide a decomposition of the total error on the flux densities into 
three parts: 

- Part of the error is the statistical error on the flux determina¬ 
tion, independent for the different wavelengths. 

- Part of the error is systematic and due to the flux difference 
between the two nod spectra and is correlated for all of the 
wavelengths over the SF range on one side and the FL range 
on the other side. 


The third component of the error in IRS spectra is the cali¬ 
bration error, and this error is correlated between the two SF 
and LF modules since the same calibrator star was used for 
both modules. According to [Lebouteiller et ak ( |201 l| l, the 
global IRS calibration is better than the 2% level. 


MIPS: For the DGS MIPS photometry, Bendo et al. (2012 1 
used a 4% calibration error at 24 pm (|Engelbracht et aL|2007] l, 
10% at 70 pm ^Gordon et al. |2007| l and 12% at 160 pm 
( Stansberry et al.|2007 1 . For KINGFISH, the same uncertain¬ 
ties were used by [Dale et al.| P007| l except for MIPS 70 pm 
where they adopted a 7% calibration error. According to the 
MIPS Instrument Handbook, the calibration of the MIPS 160 
pm band has been done using the 24 pm and 70 pm observations 
of asteroids. We can thus consider that the calibration errors for 
MIPS 24 pm and MIPS 70 pm are both correlated with MIPS 
160 pm. 


WISE; WISE calibration has been performed on stars and is 
tied to Spitzer calibration according to Jarrett et'ar| ( 201 1) . For 
each wavelength, the correlations between the bands can be sum¬ 
marised this way: 


- The WISE 3.4 pm calibration error is decomposed in an in¬ 
dependent part, proper to WISE, of 2.4 %, and is correlated 
with the IRAC 3.6 pm band. 

- Similarly, the WISE 4.6 pm error has an independent part of 
2.8% and is correlated with the IRAC 4.5 pm band. 

- The WISE 12 pm error has an independent part of 4.5% and 
is correlated with the IRS-SF/FF modules. 

- And finally, the WISE 22 pm error has an independent part 
of 5.7% and is correlated with the MIPS 24 pm band. 


IRAS; According to the IRAS explanatory supplemen f^ the 
calibration of IRAS has been tied to the |Rieke et al.| ( |1984| l 
ground-based photometric system at 12 pm. The three IRAS 
bands at 12, 25 and 60 pm have been calibrated using stellar 
models, and the IRAS 100 pm calibration used asteroids. The 
relative uncertainties relative to the ground-based 12 pm are 2%, 
5% and 5% for IRAS 12, 25 and 60 pm respectively, indepen¬ 
dent from band to band. The absolute uncertainty on the 12 pm 
flux density is 4%, in common and correlated between the three 
bands. The uncertainty at 100 pm is 10%, and is not correlated 
with any of the other IRAS bands. 


3.4. Comparison with a singie modified biackbody 


We conducted a first study of the dust properties in the DGS 
using a modBB model to describe the dust emission in |Remy 


Ruyer et al. (20131. Although very popular in the literature espe¬ 


cially when limited data is available, a modBB model assumes 
a single temperature for the dust grains, and this affects the re¬ 
sulting dust properties. Fig. shows the ratios between the dust 
masses estimated with a modBB model, Mbb, and with our semi- 
empirical SED model, for both DGS and KINGFISH sam¬ 
ples, as a function of metallicity. 

For this comparison, the modBB masses have been updated 
from Remy-Ruyer et al. (20131 to be consistent with the updated 
Herschel data. We also fixed the emissivity index (3 -2.0 to be 
consistent with the effective emissivity index of our dust model. 


This document is available at: http://lambda.gsfc.nasa.gov/ 
product/iras/docs/exp.sup/. 
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12 + log(0/H) 


0.78 0.98 1.00 1.03 1.04 1.06 1.10 1.21 1.32 1.56 2.07 

'^BB/'^dust 


Fig. 5. Ratios of the dust masses estimated with a modBB model and with a semi-empirical SED model, for the DGS 

(crosses) and KINGFISH (downward triangles) samples, as a function of the metallicity. The colour codes the ratio between the 


dust temperatures estimated from modBB fits and the SED model (assuming Tmw - 19.7 K, Planck Collaboration et al.|2014 
text for details) rBB/Tdust- 


We otherwise use the same wavelength range and fitting proce¬ 
dure as in| Remy-Ruyer et S] ( |2013| l. The opacity at the reference 
wavelength, /f(/lo) = 4.5 m^kg ‘ at 100 /rm, used to normalise the 
modBB model is derived from the opacities of the dust compo¬ 
sition adopted in the model ( Galliano et al.|20TT i. We find that 
the modBB model almost systematically underestimates the dust 
mass compared to a semi-empirical SED model, with a median 
ratio of 0.48 but also ratios sometimes as low as 0.1 (Fig. 1^, and 
that this underestimation does not depend on the metallicity of 
the source. 


We also compare the temperatures derived from a modBB 
fit, Tbb, to the average dust temperatures derived from the SED 
fits, rduist- We can estimate Tjust directly by integrating over 
T — TmwX in Eq[T] with^-2.0 for our m odel. We adopt 

Tmw= 19.7 K from [Planck Collaboration et al.| ( |2()14| l. We see 
from the colour coding on Fig. that the more T^b is overes¬ 
timated compared to Tdust, the more the modBB underestimates 
the dust mass compared to a more complex dust modelling. 

When modelling a galaxy with a broad distribution of equi¬ 
librium temperatures, the modBB model will tend to average 
over all temperatures to best fit the peak of the SED. The modBB 
will get biased towards higher dust temperatures, and thus lower 
dust masses, as seen in Fig. This also applies, to a smaller ex¬ 
tent, in the case where the dust emission is dominated by dust 
heated by the diffuse ISRF because the two big grain popula¬ 


tions (graphite and silicate) do not reach the same equilibrium 
temperatures. 


Other studies have compared dust masses from modBB to 
those from more complex dust modelling using the jPraine & Lij 
(2007 1 dust SED model (hereafter DL07), and report an under¬ 
estimation of 10 - 20% ( Magrini et al.|20lT} Bianchi|2013 1 . The 
grain mixture from [Zubko et al.j ( |20()4| l used in our modelling 
have a slightly different emissivity compared to the DL07 model 
(i.e., size distribution effect), causing the difference with the re¬ 
sults obtained here. 


3.5. Influence of the ISRF description 


We describe the starlight intensity distribution by a power-law 
component only ( Dale et al.[2001| l and not with an additional dif¬ 
fuse component as is sometimes done by other studies ( [Draine| 
[et al.|2007|[Aniano et al.|2012[[Dale et al.|20T2) l. These studies 
were dedicated to the modelling of metal-rich spirals. Our study 
is motivated by the modelling of the low-metallicity sources and 
adding a diffuse component to the starlight intensity distribu¬ 
tion would not be appropriate. Moreover, Galliano et al. ( |201 Ij l 
showed on the low-metallicity EMC that the additional diffuse 
component was not necessary. Additionally, this diffuse compo¬ 
nent would add yet another free parameter, which would be diffi¬ 
cult to constrain for the few dwarf galaxies not detected beyond 
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160 /tm. We apply this starlight intensity parametrisation to the 
KINGFIHSH galaxies as well, to have a consistent and homoge¬ 
neous modelling approach for the whole sample. Nonetheless, 
we check that the two descriptions for the starlight intensity dis¬ 
tribution (power-law only or power-law -H diffuse) give the same 
ALdu.st and {U) for the KINGFISH galaxies. We have a median ra¬ 
tio of Mdust(power-law-i-diffuse)/Mdust(power-law only) = 1.03 + 
0.15, and (t/)(power-law-i-diffuse)/(t/)(power-law only) = 1.00 
+ 0.12. Given that the average error is 20% on Mdust and 21% 
on {U), we have a good agreement between the two parametri- 
sations of the starlight intensity. 

Note that we use the same ISRF for all of the galaxies for 
consistency, not taking into account the fact that the ISRF is ex¬ 
pected to be harder in low-metallicity galaxies (Galliano et al. 
2003 2005 Madden et al.|2006l l. At hxed energy density, i.e. U, 


and given that the absorption efficiency, Qabs, is not a strongly 
varying function of wavelength in the UV/visible regime, the 
hardness of the radiation field would not greatly impact the emis¬ 
sion from the big grains in thermal equilibrium. Under these as¬ 
sumptions, increasing the hardness of the ISRF in dwarf galax¬ 
ies would only affect the emission spectrum of the stochastically 
heated grains. 


3.6. Influence of the dust composition: amorphous carbons 

Several studies suggested that amorphous carbon grains pro¬ 
vide a better description of carbonaceous dust in the ISM than 
graphite. For example, Serra Dlaz-Cano & Jones P008 1 showed 
that the erosion of amorphous carbon grains in shocks is more 
efficient than for graphite grains and matches better the observa¬ 
tions of high fractional abundances of carbon in the gas phase of 
shocked regions. They conclude that even if the carbon dust is 
in form of graphite when injected in the ISM, it is unlikely that 
it will remain graphitic as it evolves in the ISM and is subject to 
erosion or ion irradiation. 

We compute another set of dust masses for our sample using 
amorphous carbons instead of graphite grains in our full SED 
model, and keeping the same carbon mass budget in the grains. 
The amorphous carbons optical properties are taken from Zubko 
et al. (19961 (see Galliano et al. 2011|l. We use exactly the same 


procedure and options (e.g., additional MIR modBB or not) to be 
able to directly compare the graphite (noted with [c^]) and amor¬ 
phous carbon (noted with [Ac]) dust masses (reported in Table 

0 . 

We hnd that the amorphous carbon dust masses are about 
2.5 times lower than the graphite dust masses. This is because 
the amorphous carbon dust is more emissive in the submm do¬ 
main and needs less dust mass to account for the same lumi¬ 
nosity. Galliano et al.|(20TT]l showed that this dust composition 


was better suited for the low-metallicity EMC. However, we do 
not hnd any dependence of Mdust [Ac]/Afdust [Cf] on metallicity, nor 
with any other galaxy property. This is because we hxed the dust 
composition a priori in our models. Regarding the other param¬ 
eters of the ht, Ltir, ctU, and /pah do not vary signihcantly be¬ 
tween the two dust compositions, while {U)[Ac]/{U)[Cr] ~ 2.0. At 
hxed SED shape, i.e. hxed average equilibrium temperature, the 
radiation held intensity required to heat the amorphous carbon 
grains is higher as they have a higher emissivity in the submm 
and absorb more in the UV/visible regime. 

In most studies similar to ours where dust SED modelling is 


used to derive and interpret dust properties (e.g.. 

Jraine et al. 

2007 

Galliano et al.l 20081 IGalametz et aL||2011 

|Dale et al. 

2012 

Cortese et al. l2012a| Ciesla et al. 2014| | 

lemy-Ruyer 

et al. 

2014a[l the models mostly use graphite for the carbona- 


ceous component. So to ease the comparison with these various 
studies, we use Mdust[Gr] in the rest of this paper. As the ratio 
Mdust [Ac]/Mdust [Gc] is mostly independent of any of the galaxy 
properties tested here, using Mdust[Ac] instead would systemati¬ 
cally shift the trends presented in the following sections without 
affecting the general conclusions. 


3.7. Influence of the wavelength coverage and submm 
excess 


The importance of submm observations longwards of 160 pm 
has been shown by |Gordon et al. ( |2010|l ; [Galametz et al.| ( |20lT| l; 
[Dale et ah] ( |2012 l; ^iesla et ah ( |2014| l. We find fhaf the dust 
masses estimated without Herschel constraints are overesti¬ 
mated for high-metallicity galaxies (i.e., 12-i-log(0/H) > 8.2) 
and underestimated for lower metallicity galaxies in agreement 
with the conclusions of the previously mentioned studies. 


Excess emission at submm wavelengths (850 - 870 ;um), ap¬ 
pearing around 400 - 500 pm, and presently unaccounted for by 
the standard dust emission models, has been observed in some 


low-metallicity systems or low-mass spirals (e.g.,|Galliano et al. 

2003 

2005 

IDumke et al.| 20041 

Bendo et al. |20061 IGalametz 

et al. 

2009r 

20121 IZhu et al.|2009 

Bot et al. 120101 IGrossi et al. 

2010 

Dale et al. |2012l Remy-Ruyer et al. 

120131 ICiesla et al. 

2014 

Galametz et al.||20141 IGordon et al.| 

20141 IGrossi et al. 

2015 

. Several hypotheses have been made to explain this excess 


emission, but so far remain inconclusive. As the submm excess 
is not the main focus of this paper, we do not discuss them here. 
The presence of this submm excess can impact the dust mass, as 
a large amount of cold dust could be contributing to this submm 
emission. 


We adopt the same dehnition of the excess as that in|Remy- 


Ruyer et al. (2013|l, i.e., an excess is present if the residual at 


500 pm is greater than its error bar (Icr criteria). Note that the 
model used in |Remy-Ruyer et^ (2013 1 , against which the 500 
pm emission is compared, is a modBB with a hxed /? of 2.0, and 
that the SPIRE Hux densities have been updated to match the 
latest SPIRE calibration and beam areas for this study (see Sect 
2.2 1 . The hndings of Remy-Ruyer et al. (2013|l (identihcation of 


the galaxies with excess and number of excess sources) are thus 
not directly comparable to our results. 


Eight galaxies exhibit an excess at 500 pm: Haro 11, 
HS 0052-H2536, Mrk 930, NGC 1569, NGC 625, NGC 337, 
NGC 3049 and NGC 4631. Eor these eight excess galaxies (ex¬ 
cept one, NGC3049) including or omitting the 500 pm point 
in the ht yields dust mass estimates that are consistent within 
the errors, because the 500 pm excess is quite small (< 1.5 cr). 
However, omitting the 500 pm point in the ht for all of the sam¬ 
ple leads to underestimated dust masses (by a factor 2 to 4) for 
15% of the 78 galaxies detected at 500 pm, because part of the 
coldest component of the dust is not fully taken into account. 


Eor the purpose of our study, we thus need to include the 500 
pm point in the model to get an accurate estimation of the dust 
mass, and this estimation will not be biased for the eight sources 
with this small submm excess at 500 pm. 


4. Dust properties 

In this section, we scrutinise the dust properties derived from our 
realistic dust modelling over the whole IR wavelength range, and 
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present results for Mdus 0 Ltir, /pah, (U) and all. The median 
error on each parameter together with the spanned range of val¬ 
ues and the median values are shown in Table[T]for both samples. 

For each parameter we examine the relation with metallicity 
expressed as 12H-log(0/H), the star-formation activity traced by 
the sS FR, estimated in Sect. |2. 6. 2] and the stellar mass Mstar- For 
simplicity in the notations and to ease the reading, we use Z for 
12-i-log(0/H), even if Z is not strictly equivalent to the oxygen 
abundance. 

We compute Spearman rank correlation coefficient^for ev¬ 
ery pair of parameters and present them in Table ^to determine 
the most significant relations to investigate. For the 98 galaxies, 
we determine that the correlation (or anti-correlation) is signifi¬ 
cant if Ipl > 0.37 at a significance level of 0.01% (i.e., the proba¬ 
bility that the two variables are motonically correlated is greater 
than 99.99%). 


Table 2. Dust properties; correlations from Spearman rank coef¬ 
ficients. 


Param 

12+log(0/H) 

sSFR 

-^star 

12+log(0/H) 

1.00 

-0.78 

0.84 

^dust 

0.78 

-0.72 

0.92 

Ltir 

0.64 

-0.44 

0.87 

/pah 

0.67 

-0.77 

0.70 

(U) 

-0.47 

0.66 

-0.35 

o-U 

-0.52 

0.71 

-0.44 

^dmt/^star 

0.02 

-0.04 

-0.06 


4.1. Dust mass 


The dust masses in our sample cover a range of more than five 
orders of magnitude. The median value is 4.4 x 10^ Mq for 
the DGS sample and 1.9 x 10^ Mq for the KINGFISH sample, 
about two orders of magnitude above the median dust mass for 


the DGS, in agreement with the findings of Remy-Ruyer et al. 
( |20T3l l. 

The median error on the dust mass is ~20% (Table [T]). This 
error only takes into account the uncertainties in the observa¬ 
tions. No systematic modelling uncertainty is included due to the 
difficulty of determining it. However, with the different tests per¬ 


formed in Sects. 3.4 3.6 and 3.7 where we changed the model. 


the dust composition or the submm wavelength coverage, we 
can say that we have a conservative systematic modelling uncer¬ 
tainty of a factor 2 to 3 on the dust mass estimate. 

As seen in Table the stellar mass gives the strongest cor¬ 
relation with the dust mass (p(Mdust, 4Lstai)=0.92), and not the 
metallicity: p(Mdust, Z) = 0.78. The correlation between stellar 
mass and metallicity is also stronger than between dust mass 
and metallicity, p(Mstar, Z) = 0.84. This holds whether we con¬ 
sider the full sample, low-metallicity sources only or metal-rich 
sources only. This simply translates a scaling effect where a 



|et al. 12014) 1. The correlation between the dust mass and metallic¬ 
ity is thus a direct consequence of the mass-metallicity relation 
). This is illustrated in Figure]^ The best fit 

relation gives: 

log(Mdust) = (-15.0 + 0.9) + (2.6 + 0.1) x (12 + log(0/H)) (3) 
with a dispersion of 0.74 dex from the relation. 


(Tremonti et al. 2004 



Fig. 6. Dust masses (top), TIR luminosities (bottom) for the DGS 
(crosses) and KINGFISH (downward triangles) samples as a 
function of metallicity, colour coded with Mstar- The best power- 
law fit is indicated as a black line, together with the Icr disper¬ 
sion as dotted lines for the two panels. The distribution of each 
parameter is indicated on the side of each panel for both samples: 
solid line for DGS and dashed line for KINGFISH. 


4.2. Total infrared luminosity 
4.2.1. Estimating Ltir 


Several definitions exist to derive Ltir. We chose to estimate 
Ltir by integrating the best-fit modelled dust SED, L 7 -/«[S£D], 
(i.e., not including the stellar continuum) between 1 and 1000 
pm. Some studies use [8-1100] pm or [3-1100] pm intervals for 
the integration, but using any of the two other intervals gives 
equivalent results (see Table |^. 

Comparing Lrir and Ltir." We also compare the FIR lumi¬ 
nosities, LpiR over the interval [50-650] pm used in Remy-Ruyer 


|et al.| ( |20T3] l, with Ltir. We find that the Lrir accounts for about 
65+13% of the TIR luminosity, and gets lower as the peak of 
the SED broadens ((p(Lpir/Ltir, crU) - -0.83)). As the peak 
of the dust SED shifts to shorter wavelengths and broadens, 
the [50-650] pm EIR interval does not capture the bulk of the 
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Table 1. DGS and KINGFISH SED parameters. 


Parameter 

Median err 
[%] 

min 

DGS 

median 

max 

min 

KINGFISH 

median 

max 

12±log(0/H) 

26* 

7.14 

IZwl8 

7.93 

8.43 

He2-10 

7.54 

DD0154 

8.35 

8.77 

NGC1316 

Log(Msbu.) [Log(Mo)] 

47 

6.51 

HS0822±3542 

8.58 

10.62 

UM448 

7.24 

DDO053 

10.44 

11.68 

NGC1316 

Log(SFR) [Log(Moyr-')] 

22 

-2.21 

UGC4483 

-0.66 

1.40 

Haro 11 

-2.36 

M81dwB 

-0.38 

0.90 

NGC2146 

l.og(sSFR) [Log(yr-')] 

52 

-11.00 

Tol0618-402 

-9.16 

-8.13 

SBS0335-052 

-12.08 

NGC1316 

-10.35 

-9.35 

DDO053 

Log(Mdi,s,) [Log(MQ)] 

19 

2.44 

UGC4483 

5.69 

7.70 

UM311 

4.04 

DDO053 

7.28 

8.16 

NGC5457 

Log(iTiR) [LoglLe)] 

3 

6.42 

UGC4483 

8.98 

11.30 

Haro 11 

6.65 

M81dwB 

9.81 

11.13 

NGC2146 

/pah[/paho] 

17 

0 ** 

0.21 

0.66 

NGC4449 

0 ** 

0.90 

1.64 

NGC4254 

Log«(/» [Log(Uo)] 

21 

-0.51 

NGC6822 

1.30 

3.44 

SBS0335-052 

-0.62 

NGC4236 

0.29 

1.41 

NGC1377 

Log(o-£/) [Log(Uo)] 

40 

0.17 

UM311 

2.50 

4.77 

SBS0335-052 

-0.59 

NGC2841 

0.58 

3.14 

NGC2798 

Log(Md„s,/Ms,ar) 

51 

-4.53 

SBS0335-052 

-3.08 

-1.79 

Pox186 

-4.02 

NGC5866 

-2.96 

-2.02 

NGC4236 


*: corresponding to 0.1 dex, see |Madden et al. ( |2013^ . 

**: 0 for the minimum value means that PAHs are not detected in at least one galaxy in the sample. 
Note: /pAHo = 4.57 % in our model and U© = 2.2 x 10“^ W.m“^. 


emitted dust luminosity anymore. For very broad SEDs, we can 
miss up to 70% of the TIR luminosity. The most extreme case is 
SBS 0335-052 where Lfir/Etir= 10%, consequence of its very 
peculiar SED peaking around 30 /tm. 


Table 3. Ltir comparisons. 


X 


Median ratio 


Definitions 
TIR: Ltir[8-iioo] 
TIR: Ltir[3-iioo] 
FIR: LpiR[50-6501 


X / Ltir[1-1000] 
I.00±0.02 
1.05 ±0.08 
0.65 ±0.13 


Galametz et al. (|2013b 

7D - 


x/L- 


TIR[3-1100] 


1.13 ± 0.29 


100 

24-160 
24 - 70 - 160 
24 - 70 - 100 - 160 


0.72 ± 0.27 
1.01 ±0.35 
1.01 ±0.14 
1.00 ±0.13 


Madden et al. 1 2013 i;|Kennicutt et al. j201l[ x / Ltir[24-7o-i60.g7^ 


^-^lR[pve-Herschel] 


1.13 ±0.29 


Prescriptions for the DGS: Galametz et al. ( |2013b) l explored 
various Ltir calibrations from Spitzer and Herschel bands, using 
the KINGFISH sample. We test all of the possible relations using 
24, 70, 100 and 160 pm, and look at the ratio of the estimated 
Ltir over the Ltir derived from our best-fit SEEp^to determine 
which relations are the most appropriate for the DGS sample. 

We find that for the DGS galaxies, the best monochromatic 
relation is using 70 pm with a median ratio of 1.13 + 0.29 (Table 
[^. Using 70 pm tends to overestimate the Ltir for the DGS, 


while using 100 pm tends to underestimate the Ltir (Table [^. 
This is due to the higher Lvo/Lioo colour in dwarf galaxies than 
in more-metal rich galaxies over which these relations have been 
calibrated. The best relation using two bands is with 24 and 160 
pm. Increasing the number of bands gives similar results, with 
decreasing scatter. The smallest dispersion using three bands is 
achieved by combining 24, 70 and 160 pm. Thus to estimate 
Ltir for galaxies with similar metallicity and star-formation ac¬ 
tivity, we recommend using the 70 pm, the 24-160 pm, or the 
24-70-160 pm calibration depending on the number of available 
constraints. 

Pm-Herschel estimates of the Ltir derived with Spitzer 
fluxes and the [Dale & Helouj (|2002|) formula (LTiRipre /ferac/ic/i) 
from |Madden et al.| ( |2013p Wd jKennicutt et d^p011] l are con¬ 
sistent with the Ltir derived here. All of these Ltir comparisons 
are summarised in Table and show that the TIR luminosity is 
quite a robust parameter. 


4.2.2. Ltir and metallicity 

Ltir is presented in Fig. [^(bottom panel) as a function of metal¬ 
licity. The Ltir in both samples cover five orders of magnitude, 
and the low-metallicity dwarf galaxies are about 6.5 times less 
luminous in the IR than the more metal-rich environments. As 
for the dust mass, Ltir is strongly correlated with the stellar 
mass,p(LTiR, Mstar)=0.87, due to scaling effects. The correlation 
between Ltir and metallicity (p(Ltir, Z)=0.64) is also a conse¬ 
quence of the mass-metallicity relation. As seen in Fig. the 
relation of Ltir with metallicity is more dispersed than between 
the dust mass and the metallicity. The best fit relation gives: 

log(LTiR) = (-6.6 ± 0.2) ± (2.0 ± 0.1) x (12 + log(0/H)) (4) 


** For this comparison we use L^ir integrated over [3-1100] pm to 
match the interval used in ( 


Galametz et al. 


1 2013b I. 


with a dispersion of 0.79 dex from the relation. From the con¬ 
struction of the model, Ltir is directly linked to the product 
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MdustX{U) at first order. And while Mdust is correlated with 
metallicity, (U) is anti-correlated with metallicity, which results 
in a flattened correlation with larger scatter between Ltir and 
metallicity. Ltir is not correlated with (U) Ltir) = - 

0.05), indicating that it is the amount of dust that primarily drives 
the luminosity. 


4.3. PAHs 

Lower PAH abundances in low-metallicity systems than in more 


metal-rich galaxies have previously been observed (e.g. Madden 

2000 

Boselli et al. 2004| Engelbracht et al. |2005| Wu et al. 

2006 

O’Halloran et al. 

[20061 [Draine et al. 2007[ Smith et al. 

2007t 

Gordon et al.|2008l Galliano et al.|2008| Wu et al.|2011| 

Ciesla et al. 

2014 

1. Erom our model we obtain an estimate of 


the mass fraction of PAHs for our sources, including the 19 low- 
metallicity galaxies with new IRAC and IRS data. We have an 
average error of 17% on /pah, and a very good agreement be¬ 
tween our /pah to the PAH mass fraction derived by Dale et al. 
( 2012| l for the KINGFISH sample using the Draine & Li 1 2007| ) 
model (median ratio of 1.03 + 0.24). For the DGS galaxies with 
featureless IRS continuum, we fit the model to the observations 
again, leaving /pah free, and use this value as an upper limit, 
reported in Table [9] 

In our sample, the mass fraction of PAHs covers two orders 
of magnitude, and are detected over 1.1 dex in metallicity down 
to ~0.08 Zq, expanding the range explored by previous works. 
From Table /pah correlates very well with sS FR and metal¬ 
licity (p(/pAH, sSFR) - -0.77, and p(/pah, Z) = 0.67, not includ¬ 
ing the upper limits), with a PAH mass fraction decreasing with 
decreasing metallicity and increasing sS FR, confirming results 
from the previously mentioned studies. This is shown in Fig. 
The best-fit relations of /pah with metallicity and sS FR give; 


/ log(/pAH) = (-11.0±0.3) + (1.30±0.04)x(12-Hog(0/H)) 
1 log(/pAH) = (-5.5 ±0.1)-(0.53 ± 0.01) xlog(sSFR) 


(5) 


with a dispersion of 0.35 dex around the relation with metallicity, 
and of 0.38 dex around the relation with sS FR. 

The lower abundance of PAHs in dwarf galaxies is not con¬ 
trolled by a single parameter but rather arises from a joint effect 
of both low metallicity and high sS FR. The higher sS FR results 
in a harder and more intense galaxy-wide ISRF. Combined with 
the lower dust attenuation, PAHs are efficiently destroyed by 


hard UV photons but also shocks and cosmic rays (e.g.. Madden 
|et al.|2006{|Engelbracht et al.|2008||Micelotta et al.|2011| l. Due 

to the small physical size of dwarf galaxies, PAHs are also sub¬ 
ject to destruction by SN shock waves on galaxy wide scales 
(|0’Halloran et al.|2006||Micelotta et al.|2010 1 . In addition, lower 


C/O ratios in dwarf galaxies (Garnett et al. |1995| |1999j l mean 
that less material is available in the ISM to form the PAHs than 
to form the oxygen-rich silicate grains which make the bulk of 
the dust mass. The deficiency of PAHs in low-metallicity galax¬ 
ies can also be explained by the delayed injection of carbon dust 
into the ISM by AGB stars ( [Galliano et al.|2008] l. 


4.4. Temperature distribution 


Remy-Ruyer et al^ (|2013|l showed that dwarf galaxies harbour 


warmer dust and present a potentially broad dust SED peak. 
Warmer dust in dwarf galaxies had been discovered first with 
IRAS (e.g. |Helou|[T98^ [Hunter et al.|[T989l [M^se & Israelj 
[1994[ [Galliano et al.[[2003[l, and this was confirmed later with 
Spitzer (e.g. Galliano et al. [2005] [Rosenberg et al.[2006[[Cannon[ 


[et al.[200^[Galametz et al.[2009| l. More recent studies have con¬ 
firmed that low-mass galaxies have broader IR SED peaks (e.g. 
[Boselli et al.[[2012l [Smith et al.[[2012al [Ciesla et al.[[2014) l, or 
show a flattening of their EIR SED slope compared to more mas¬ 
sive systems ( [Boselli et al. 2010t Cortese et al.[ 2014[ l. These 
studies also showed that the most actively star forming galaxies 
were the ones with the warmest dust and the broadest SED peak 
( [Boselli et al.|2010[[a^a et al.|2014[ ). 



1 10 100 1000 
Wavelength [/zm] 


Fig. 8. Example of different SED shapes, with increasing broad¬ 
ness around the EIR peak for increasing crU. NGC 0628, 
NGC 1140 and SBS 1415±437 have log(o-t//Uo)= 0.4, 2.1 and 
3.0, respectively. The location of these three galaxies is indicated 
inEig.g 


A temperature distribution is necessary to properly describe 
the dust emission (Sect. [3^ , and this dust temperature distribu¬ 
tion is directly linked to the range of starlight intensities. At/, 
to which the dust is exposed. However, to avoid being affected 
by potential degeneracies between the starlight intensity distri¬ 
bution parameters (At/, t/min and a) we consider the standard 
deviation of the starlight intensity distribution, crU, to describe 
the range of starlight intensities to which the dust is exposed. 
Fig. [^illustrates how crU traces the width of the SED peak on 
three examples. Eig. [9| shows the range of starlight intensities 
to which the dust is exposed, crU, as a function of the average 
starlight intensity {U}. 

As explained in Sect. [3.2[ we added a MIR modBB to get 
a better match to the observed MIR SED for eleven galax¬ 
ies. We add a delta-function describing this additional single 
temperature component to the U distribution in Eqs 0 and [^ 
6(U - I/mir.bb), with I/mir.bb - (Tmir,bb/Tmw)^- ( U ) covers 
a range of four orders of magnitude equivalent to ~ 50 K in dust 
temperature, from 12 K to 64 K in the DGS; and from 12 K to 
30 K in the KINGEISH sample. The median T^ust is 26 K for the 
DGS and 20 K for the KINGEISH sample, consistent with low- 
metallicity galaxies harbouring warmer dust than more metal- 
rich environments ([Remy-Ruyer et al.|2013|l. The median crU is 
290 Uo for the DGS and 4 Uq for the KINGEISH sample, re¬ 
flecting the broader SED peak for low-metallicity environments 
compared to more metal-rich sources. 

Both {U) and crU strongly correlates with sSFR (Table 
1^, although the correlation is stronger with crU. The correla¬ 
tion of both parameters with metallicity is weaker than with 
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Fig. 7. PAH mass fractions for the DGS (purple crosses) and KINGFISH (orange downward triangles) samples as a function of 
metallicity (left panel) and sSFR (right panel), /pah is expressed in units of /pahg^ with /paho = 4.57%. The distribution of /pah 
is indicated on the side for both samples: plain purple line for DGS galaxies and dashed orange line for the KINGFISH sample. In 
each panel, the best power-law fit is indicated as a solid line, together with the Icr dispersion as dotted lines. 
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Fig. 9. crll as a function of the average starlight intensity (U) for the DGS (crosses) and KINGFISH samples. For the KINGFISH 
sample, we distinguish the sources that have been classified as “star-forming” (SF, downward triangles) or as “AGN” (stars) by 
Kennicutt et al. ( 201 1[ ). The average dust equilibrium temperature is shown on the top axis. The colour codes the sSFR of the 
sources, in units of Gyr'F We indicate the location of the two most metal-poor galaxies, I Zw 18 and SBS 0335-052, and of the 
three galaxies of Fig. 
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sSFR (p(Z, crU)--0.52, p(Z, {U})--0.47), and is a side ef¬ 
fect of the strong correlation between sSFR and metallicity 
(p(Z, sSFR) = -0.78). This holds whether we consider the full 
sample or low-metallicity sources on one side and metal-rich 
sources on the other side. The sS FR thus seems to be the pa¬ 
rameter driving the dust SED shape. 

In Fig. ^ we clearly distinguish two clusters of points, one in 
the lower eft corner of the plot, mostly metal-rich KINGFISH 
galaxies, with low (U) and low crU (“cold and narrow” SEDs), 
and one mainly composed of low-metallicity galaxies, with 
higher (U) and crU (“hot and broad” SEDs). However some 
high-metallicity sources can be found in the “hot and broad” 
group, and vice versa. This is due to their higher (or lower) 
sS FR as shown by the colours on Fig. For example, the spi¬ 
ral galaxy NGC 0628 has the same metallicity than the star¬ 
forming dwarf galaxy NGC 1140 (Z = 0.5 Zq) but very different 
SED shapes (Fig.[^ due to their different star-formation activity: 
logisSFR/Gyr^) = -1.4 for NGC 0628 and -0.8 for NGC 1140. 
This confirms that the sS FR is the parameter determining the 
dust SED shape, and that metallicity only plays a secondary role 
as also noted by da Cunha et al. ( |2008| l; Smith et al. ( 2()12a| l. 
Low-metallicity sources extend the trend outlined by metal-rich 
galaxies to warmer temperatures and broader temperature distri¬ 
butions due to their higher star-formation activity in a smooth 
transition rather than in a sharp change. This has been noted by 
Cortese et al. ( 2014| l on a smaller range of metallicities (0.6 dex). 

In more active galaxies, the dust spans a wider temperature 
range translating into a broader dust SED peak. In a galaxy 
undergoing an active phase of star formation (e.g., with high 
sS FR), the ISM will be clumpier as a large number of embedded 
star-forming clumps are spread all over the galaxy. This clumpier 
ISM allows for a wider equilibrium temperature distribution of 
the dust grains, skewed towards higher dust temperatures (and 
thus higher crU and (t/)). Evidence for this clumpier structure 
of the ISM in dwarf galaxies can be directly seen from the re¬ 
solved UV-to-mm observations of the EMC and SMC, and has 
also been suggested by Cormier et al.| ( 2015| l from a detailed 
study of the DGS FIR hne-stracture cooling lines. This irregular 
ISM stracture due to feedback processes related to star formation 
and supernovae events result in an ISM in which the dust grains 
are exposed to a range of stellar populations and thus show a 
larger distribution of dust temperatures. The dwarf galaxies in 
our sample contain warmer dust primarily because they are in¬ 
tensively forming stars. Metallicity has a secondary impact on 
the dust temperature, as the low dust attenuation allows to heat 
the dust deeper within the molecular clouds. 

In Fig. 1^ we have also split the KINGFISH sample in 
the “SF” and “low-luminosity AGN” groups. The warmest 
KINGFISH galaxies (i.e., with high (t/)) are mostly the 
KINGFISH low-luminosity AGNs. It might be possible that, de¬ 
spite a small contribution to the total luminosity, the emission 
from the central AGN is powerful enough to impact the dust 
heating on global galaxy scales and to increase the average dust 
temperature, thus to have an impact on the global shape of the 
dust SED. Warmer dust in the presence of an AGN had already 
been seen in Kirkpatrick et al. (2012|l. 


The most famous picture of galaxy evolution is the widely 
studied mass-metallicity relation. The ISM matter lifecycle im¬ 
plies that the metallicity of a galaxy increases as the galaxy 
evolves through several cycles of star formation, building up its 
stellar mass ([Lequeux et al.||1979[ Tremonti et al.||2004 Kirby 


eTaLjlTOlA) [Cook et al.||2014| [Gongalves et al.||2014| l. We 


saw 

that the mass-metallicity relation drives all of the extensive quan¬ 
tities, Mdust and Ltir, while the intensive quantities /pah, {U) 
and crU, are mostly driven by sS FR. This is shown schemati¬ 
cally in Fig. From this hgure, is is clear that metallicity only 
plays a secondary role in shaping the dust properties. Remember 
however, that our sample is dominated at high metallicities by 
spiral galaxies and by star-forming, gas-rich dwarf galaxies at 
low metallicities. The correlations presented here thus may suf¬ 
fer from some selection bias. The dust-to-stellar mass ratio being 
a normalised (and thus intensive) quantity enables us to look at 
the dust mass build up with respect to the stellar mass. Fig. 11 
shows the dust-to-stellar mass ratios for our sample as a function 
of sSFR. 

The dust-to-stellar mass ratio has been extensively stud- 
ied previously (e.g. |da Cunha et al.||2010[ [Skibba et al.||2011t 


[Bourne et al.||201^ Cortese et al.||2012bHSmith et al.||2012c] l. 

These studies have shown that the dust-to-stellar mass ratio de¬ 
creases for increasing stellar mass (or metallicity) and decreas¬ 
ing sS FR. We do not find any correlation for the dust-to-stellar 
mass ratio with metallicity, stellar mass or sSFR. However, con¬ 
sidering each sample separately, we find a rather weak cor¬ 
relation between the dust-to-stellar mass ratio and metallicity 
(p(Mdust/M„ar, Z)=-0.44) or sSFR (p(Mdust/M„„, sSFR)^Q.41) 
for the KINGFISH sample, and still no correlations for the DGS 
sources : p(Mdust/Mstar, Z)=0.18 and p(Mdust/Afstar, sSFR)-- 
0.04. Thus the results for the metal-rich KINGFISH galaxies are 
in agreement with the findings of the previously mentioned stud¬ 
ies. However, the low-metallicity sample does not extend the ob¬ 
served behaviour of the metal-rich galaxies to higher sS FR. 

The peculiar behaviour of the dust-to-stellar mass ratios for 
the low-metallicity DGS galaxies is due to their chemical evo¬ 
lutionary stage. This can be traced by the gas-to-dust-mass ratio 
(G/D, [Remy-Ruyer et alT 2014a[ ), and is shown by the colours 
in Fig. jl 1[ Modelled evolutionary tracks from the chemical evo¬ 
lution model of [Asano et al.[p()13[ l are also shown in Fig. 
for different star-formation timescales (tsf = 0.5, 5 and 50 Gyr). 
This chemical evolution model, based on models from [Hirashita[ 
( 1999| l; [Inoue| ( [20lT] l, includes dust production by stellar sources 
(AGB stars and Type II SNe) and by dust growth processes in 


the ISM. Dust is destroyed by SN shocks. [Remy-Ruyer et ah 
( [20I4a| l showed that this model can successfully reproduce the 
observed trend between G/D and metallicity, thanks to the dust 
growth mechanism. 

At high sS FR, thus low-metallicity and low stellar masses, 
the dust content is very low giving very low Mdust/Afstar and high 
G/D. The dust production is controlled there by stellar sources 
only. Then when a critical metallicity is reached in the ISM, 
dust growth by metal accretion on the dust grains in the ISM 
becomes the major process for building up the dust mass (|Asano[ 


et al. 2013[ Zhukovska 2014 1 and the dust mass rapidly increases 


5. Towards a comprehensive view of the dust 
properties in iow-metaiiicity environments 

After studying the different relationships between the dust prop¬ 
erties and fundamental galaxy parameters, we now analyse our 
results in the context of galaxy evolution, and attempt to draw a 
consistent picture of the dust in low-metallicity environments. 


without significant consumption of the gas reservoir or star for¬ 
mation. This rapid increase of the dust mass results in an in¬ 
crease of the dust-to-stellar mass ratio and a decrease of the 
G/D ( Remy-Ruyer et al.|2014a| l. The large scatter in the dust-to- 
stellar mass ratio for the highly star forming galaxies in Fig. [TT] 
is thus due to their different evolutionary stage, as traced by their 
G/D. Then dust growth processes saturate when all the available 
metals are locked up in the grains. In the meantime, star forma- 
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Fig. 10. Schematic view of the various correlations between parameters in the DGS and KINGFISH samples. The numbers on the 
lines are the Spearman correlation coefficients reported from Table The solid lines indicates the dominating correlations between 
parameters and the dashed lines indicates several secondary correlations of interest. 


tion continues, consuming the gas reservoir and increasing the 
stellar mass. This results in a decreasing sS FR, decreasing dust- 
to-stellar mass ratio and decreasing G/D. 

The nucleus type (star-forming or AGN) also has an impact 
in shaping the dust properties (|Fritz et ar]|2006[ |Magdis et al.| 


2014 


|Ciesla et al.|2015]l, meaning that the picture drawn from 


Figs. 10 and 11 is far from being that simple. sSFR and metal- 
licity are not the only parameters involved in the evolution of the 
dust properties in galaxies. All of these fundamental parameters 
can also be affected by inflows and outflows, rendering more 
complex in practice the relations between stellar mass, star for¬ 
mation, metallicity and dust properties, and can change signih- 
cantly the position of a galaxy in the Mdust/4Lstar- sS FR plot. 

Moreover, the dust masses are estimated here with the same 
dust mixture, using silicates, graphite and PAH grains, for all 
of the galaxies. However, amorphous carbon grains could be 
used instead of graphite grains to model the carbonaceous dust 
( [Galliano et al.|201 lUGalametz et al.|2013a[ [Jones et al.||2013] 


and see Section 3.6 1 . The silicate grain properties used in our 


model are derived from empirically-normalised astrophysical 
data, and more realistic types of silicate grains could also be 
used (Kohler et al. 2014j l. Changing the dust composition results 
in different dust properties, e.g., the dust mass or dust temper¬ 
ature as we saw with amorphous carbons in Section 3.6 The 


dust composition can also vary from galaxy to galaxy, which 
would in turn impact the derived dust masses, and the observed 
trends. Thanks to the wealth of spatially resolved dust studies 
with Herschel, we know now that the dust composition and size 
distribution also vary within each galaxy ([Smith et ah [2012b 


Galametz et al.|20T^[Mattsson et al.|2014b[[Viaene et al.|2014 


Ysard et al.^m^ Dust evolves from diffuse to denser regions 

where bigger grains or aggregates can eventually form (Kohler 
et al.|2015] ), thus impacting the SED shape, and the dust param¬ 
eters. Studying the impact of a varying dust composition and 
grain size distribution within and between galaxies on the SED 
and the derived dust properties is now the next step. 


6. Conclusion 

In this paper, we provided a coherent picture of the evolution of 
the dust properties from metal-poor to metal-rich environments. 
Our sample comprises 109 galaxies, spanning almost 2 dex in 
metallicity, dominated by spiral galaxies at high metallicities 
and by star-forming gas-rich dwarf galaxies at low metallicities. 
Observed SEDs are gathered over the whole IR-to-submm wave¬ 


length range, with constraints from 2MASS, Spitzer, WISE, 
IRAS, and Herschel. The full data set is presented here for the 
DGS sample of dwarf galaxies. The dust properties (namely, dust 
mass, TIR luminosity, PAH mass fraction and dust temperature 
distribution) are derived in a systematic way using a realistic 
semi-empirical dust SED model from Galliano et al. ( 201 l| l, and 
then compared to fundamental galaxy parameters: stellar mass, 
metallicity and sSFR for a hnal sample of 98 galaxies. 

The dust mass is a critical parameter for constraining chem¬ 
ical evolution models. We showed here that different model 
assumptions could greatly impact the estimated dust mass. A 
single-temperature modBB model underestimates the dust mass 
by a factor of 2 compared to a semi-empirical SED model, even 
with careful matching of the effective optical properties in the 
modBB. The modBB model overestimates the effective average 
dust temperature thus leading to the underestimation of the dust 
mass, as also shown by Bendo et al. (2015 1 . Changing the car¬ 
bonaceous component in the dust mixture from graphite grains 
to amorphous carbon grains results in a decrease of the dust mass 
estimate by a factor 2.5. Amorphous carbon grains are more 
emissive than graphite grains, so less dust is needed to account 
for the same luminosity. 

We hnd an excess at 500 pm for eight galaxies in our sample. 
The excess is rather small (< 1.5 cr) as the transition from ther¬ 
mal dust and submm excess emission observed at longer wave¬ 
lengths occurs around 500 pm. Including the 500 ;um point dur¬ 
ing the htting procedure does not result in a drastic overestima¬ 
tion of the dust mass for the excess galaxies. However, leaving 
the 500 pm data point out of the ht results in an underestima¬ 
tion of the dust mass by a factor 2 to 4 for 15% of the galax¬ 
ies in our sample, because the cold dust component is then not 
properly constrained. Estimating dust masses is thus subject to 
non-negligible systematic modelling uncertainties. 

We present various ways of estimating Lxir in the DGS sam¬ 
ple with mono/polychromatic indicators, using the calibrations 
presented in Galametz et al. P013b| l. Eor galaxies with similar 
metallicity and star-formation activity, we recommend using the 
70, 24-160, or 24-70-160 calibrations of [Galametz et al.[P013b| l. 

The dust temperature distribution and PAH mass fraction 
are primarily driven by the sSFR, with a second order effect 
from metallicity. In our sample, PAHs are detected down to Z 
~ 1/12 Zq. Low PAH abundances in dwarf galaxies are a conse¬ 
quence of the high star-formation activity and low dust attenu¬ 
ation. This combined effect of sS FR and metallicity is also re¬ 
sponsible for the higher average dust temperatures in starburst- 
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Fig. 11. Dust-to-stellar mass ratios for the DGS (crosses) and KINGFISH (downward triangles) samples as a function of sS FR, 
colour coded by the gas-to-dust mass ratio (G/D), from |Remy-Ruyer et al.|(|2014a|l. The distribution of M^ust/AL, 

m. 


IS 


indicated on the side for both samples: solid line for DGS and dashed line for KINGFISH. Modelled evolutionary tracks from Asano 
|et al.|([20T3|l are shown for different star-formation timescales (tsf = 0.5, 5 and 50 Gyr) by the red, blue and purple lines. 


ing low-metallicity dwarf galaxies. We find a median dust tem¬ 
perature (derived from the average starlight intensity {U)) for the 
DGS sample of ~26 K, and ~20 K in KINGFISH galaxies. The 
higher star-formation activity results in a clumpier ISM allowing 
for a larger range of dust equilibrium temperatures. 


The extensive dust properties, Mdust and Ltir. are driven by 
the mass-metallicity relation, reflecting a scaling effect: the more 
massive galaxies contain more dust, are more luminous, and are 
also more metal-rich. However, the dust mass build up with re¬ 
spect to stellar content is not the same in highly star-forming, 
low-metallicity sources than in more metal rich systems. The 
dust-to-stellar mass ratios of metal-rich sources follow an in¬ 
creasing trend of Mdust/ALstar with sS FR, previously observed on 
other samples (|da Cunha et al.poTO Skibba et al.|201 1[ Cortese 


et al.|2012b|. On the other hand, for the more actively star form¬ 


ing galaxies (sSFR >0.1 Gyr“^) the trend is far less clear, with 
increasing scatter. The peculiar behaviour of the low-metallicity 
sources is driven by their chemical evolutionary stage: at low 
metallicity and high star-formation activity , the dust production 
is dominated by stellar sources only. After a critical metallicity is 
reached, the dust-to-stellar mass ratio rapidly increases because 
dust growth processes in the ISM dominate the dust production. 
Then dust growth saturates while star formation goes on, form¬ 
ing stars from the gas reservoir, resulting in a decreasing sS FR 
and Mdust/Mstar. This completes and confirms our results on the 
gas-to-dust mass ratios derived in Remy-Ruyer et al^ ( |2014a i. 
Others effects such as inflows or outflows, presence of an AGN, 
or varying dust composition between and within galaxies have 


not been considered here but can also affect our picture of the 
evolution of the dust properties in galaxies. 
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Fig. 12. DGS SEDs: The observed SEDs include the Herschel data (purple crosses) as well as any available ancillary data (in 
orange). The different symbols refer to the different instruments: Xs for 2MASS bands, stars for Spitzer IRAC and MIPS, diamonds 
for WISE and triangles for IRAS. The IRS spectrum is displayed in orange. The total modelled SED in black is the sum of the 
stellar (green) and dust (red) contributions, and eventually of the modified blackbody contribution (yellow). The modelled points in 
the different bands are the filled blue circles. 
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Fig. 12. DGS SEDs (continued). 
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Fig. 12. DGS SEDs (continued). 
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Fig. 13. KINGFISH SEDs: The colors and symbols are the same as in Fig. 
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Fig. 13. KINGFISH SEDs (continued). 






























[Lg] i/'L^ [Lg] [Lg] [Lg] uL^, [Lg] 


A. Remy-Ruyer et al.: Linking dust emission to fundamental properties in galaxies, Online Material p 7 



Fig. 13. KINGFISH SEDs (continued). 
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Table 4. DGS and KINGFISH stellar mass and star-formation parameters. 


Name 

Log[Ms,ar] 

Log[SFR] 

Logics FR] 

Log[L„„] 

Ref 


Log[Mo] 

Log[Moyr-'] 

Log[yr^'] 

Log[erg/s] 



(1) 

(2) 

(3) 

(4) 

(5) 

DGS 






Haro 11 

10.24 ± 0.21 

1.40 + 0.33 

-8.84 + 0.39 

42.45 

(1) 

Haro2 

9.56 + 0.21 

-0.09 + 0.20 

-9.66 + 0.29 

40.96 

(2) 

Haro3 

9.49 + 0.21 

-0.20 + 0.18 

-9.69 + 0.28 

40.83 

(3) 

He2-10 

9.50 + 0.21 

-0.27 + 0.14 

-9.77 + 0.25 

40.70 

(4) 

HS0017-H1055 

7.84 + 0.16 

- 

- 

- 

- 

HS0052-H2536 

10.31+0.30 

0.47 + 0.07 

-9.84 + 0.31 

- 

- 

HS0822-H3542 

6.51 + 0.17 

-1.96 + 0.26 

-8.47 + 0.31 

39.03“ 

(5) 

HS1222-H3741 

9.38 + 0.16 

0.29 + 0.16 

-9.08 + 0.22 

41.54 

(6) 

HS1236-H3937 

< 9.35 

- 

- 

40.37 

(6) 

HS1304-H3529 

8.58 + 0.31 

0.04 + 0.01 

-8.54 + 0.31 

41.27 

(6) 

HS1319+3224 

8.05 + 0.24 

-1.35 + 0.13 

-9.41+0.27 

- 

- 

HS1330-H3651 

8.95 + 0.28 

-0.21+0.03 

-9.16 + 0.28 

41.01 

(6) 

HS1442-H4250 

7.92 + 0.13 

-1.67 + 0.06 

-9.58 + 0.14 

39.56 

(5) 

HS2352-H2733 

8.70 + 0.16 

- 

- 

- 

- 

IZwl8 

7.34 + 0.21 

-1.12 + 0.09 

-8.46 + 0.23 

40.14 

(5) 

ICIO 

8.44 + 0.02 

-1.57 + 0.01 

-10.01+0.02 

39.54 

(4) 

IIZw40 

8.60 + 0.33 

0.04 + 0.09 

-8.56 + 0.34 

41.24 

(4) 

Mrkl089 

10.28 + 0.21 

0.65 + 0.16 

-9.63 + 0.26 

41.69 

(7) 

Mrkl450 

8.09 + 0.21 

-0.92 + 0.05 

-9.01+0.22 

40.29 

(5) 

Mrkl53 

9.06 + 0.21 

-1.04 + 0.14 

-10.10 + 0.25 

39.86 

(3) 

Mrk209 

7.46 + 0.21 

-1.55 + 0.13 

-9.01+0.25 

39.70 

(4) 

Mrk930 

9.44 + 0.21 

0.62 + 0.15 

-8.82 + 0.26 

41.77 

(3) 

NGC1140 

9.64 + 0.21 

-0.11+0.04 

-9.76 + 0.22 

41.04 

(3) 

NGC1569 

8.93 + 0.21 

-0.21+0.04 

-9.13 + 0.21 

41.02 

(4) 

NGC1705 

8.34 + 0.34 

-1.26 + 0.02 

-9.60 + 0.34 

39.98 

(4) 

NGC2366 

8.29 + 0.31 

-1.12 + 0.07 

-9.41+0.32 

40.10 

(4) 

NGC4214 

9.03 + 0.21 

-0.96 + 0.04 

-9.99 + 0.21 

40.18 

(4) 

NGC4449 

9.43 + 0.21 

-0.41+0.07 

-9.84 + 0.22 

40.70 

(4) 

NGC4861 

8.49 + 0.21 

-0.66 + 0.13 

-9.15 + 0.25 

40.58 

(5) 

NGC5253 

8.77 + 0.21 

-0.57 + 0.18 

-9.34 + 0.28 

40.54 

(4) 

NGC625 

8.73 + 0.21 

-1.30 + 0.06 

-10.03 + 0.22 

39.83 

(4) 

NGC6822 

8.19 + 0.23 

-1.89 + 0.25 

-10.08 + 0.34 

39.13“ 

(4) 

Pox186 

7.06 + 0.17 

-1.46 + 0.03 

-8.52 + 0.17 

39.77 

(5) 

SBS0335-052 

8.00 + 0.15 

-0.13 + 0.16 

-8.13 + 0.22 

41.08 

(8) 

SBS1159-H545 

7.84 + 0.18 

-0.90 + 0.06 

-8.75 + 0.19 

40.33 

(9) 

SBS1211+540 

7.60 + 0.21 

-1.73 + 0.29 

-9.33 + 0.36 

39.39“ 

(9) 

SBS1249+493 

8.42 + 0.39 

-0.40 + 0.11 

-8.82 + 0.40 

40.82 

(10) 

SBS1415+437 

7.80 + 0.26 

-1.20 + 0.11 

-9.00 + 0.29 

40.04 

(5) 

SBS1533+574 

9.33 + 0.21 

-0.17 + 0.13 

-9.50 + 0.25 

41.02 

(5) 

Tol0618-402 

10.43 + 0.21 

-0.57 + 0.11 

-11.00 + 0.24 

40.63 

(11) 

Toll214-277 

8.17 + 0.18 

-0.10 + 0.13 

-8.27 + 0.22 

41.14 

(11) 

UGC4483 

6.89 + 0.22 

-2.21+0.18 

-9.11+0.28 

38.61“ 

(4) 

UGCA20 

< 7.14 

-1.79 + 0.23 

> -8.93 

39.30“ 

(4) 

UM133 

< 7.99 

-1.11+0.06 

> -9.09 

40.13 

(5) 

UM311 

9.78 + 0.13 

0.38 + 0.18 

-9.40 + 0.22 

41.60 

(12) 

UM448 

10.62 + 0.21 

1.03 + 0.10 

-9.59 + 0.23 

42.07 

(3) 

UM461 

7.66 + 0.22 

-1.42 + 0.15 

-9.08 + 0.26 

39.80 

(3) 

VIIZw403 

7.21 + 0.21 

-1.89 + 0.25 

-9.10 + 0.33 

39.14“ 

(4) 

KINGFISH 






NGC0337 

10.11+0.17 

0.18 + 0.05 

-9.93 + 0.18 

41.25 

(2) 

NGC0584 

10.87 + 0.20 

-1.15 + 0.10 

-12.02 + 0.22 

- 

- 

NGC0628 

10.29 + 0.21 

-0.14 + 0.08 

-10.43 + 0.22 

40.87 

(2) 

NGC0855 

9.20 + 0.21 

-1.27 + 0.05 

-10.47 + 0.21 

39.81 

(2) 

NGC0925 

9.97 + 0.20 

-0.24 + 0.06 

-10.21+0.21 

40.89 

(2) 

NGC1097 

11.00 + 0.20 

0.48 + 0.09 

-10.53 + 0.22 

41.23 

(2) 

NGC1266 

10.18 + 0.21 

0.17 + 0.10 

-10.01+0.24 

40.19 

(3) 

NGC1291 

11.02 + 0.20 

-0.47 + 0.06 

-11.49 + 0.21 

40.64 

(4) 

NGC1316 

11.68 + 0.20 

-0.40 + 0.14 

-12.08 + 0.24 

40.29 

(13) 

NGC1377 

9.47 + 0.21 

- 

- 

- 

(3) 

NGC1404 

11.15 + 0.20 

-0.70 + 0.20 

-11.85 + 0.29 

40.59 

(13) 

IC0342 

10.91+0.04 

0.03 + 0.10 

-10.88 + 0.11 

40.73 

(4) 

NGC1482 

10.55 + 0.21 

0.46 + 0.15 

-10.10 + 0.26 

40.80 

(2) 
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Table 4. (continued) DGS and KINGFISH stellar mass and star-formation parameters. 


Name 

Log[Mstar] 

LogiMo] 

(1) 

Log[SFR] 
Log[M0 yr-‘] 
(2) 

Log[i5FR] 

Log[yr^‘] 

(3) 

Log[L„„] 

Log[erg/s] 

(4) 

Ref 

(5) 

NGC1512 

10.36 ±0.19 

-0.52 ±0.15 

-10.89 ±0.24 

40.42 

(2) 

NGC2146 

10.94 ± 0.04 

0.90 ±0.14 

-10.05 ±0.14 

41.37 

(3) 

HoII 

8.23 ±0.21 

-1.51 ±0.15 

-9.74 ± 0.25 

39.73 

(4) 

DDO053 

7.24 ± 0.32 

-2.11 ±0.08 

-9.35 ±0.33 

38.78“ 

(2) 

NGC2798 

10.36 ± 0.25 

0.39 ± 0.20 

-9.98 ±0.32 

40.95 

(2) 

NGC2841 

11.08 ±0.20 

-0.10 ±0.08 

-1L18±0.21 

40.81 

(2) 

NGC2915 

8.49 ± 0.22 

-1.78 ±0.17 

-10.28 ±0.27 

39.43 

(2) 

Hoi 

7.84 ±0.14 

-2.08 ± 0.06 

-9.92 ±0.16 

38.83“ 

(4) 

NGC2976 

9.32 ±0.21 

-1.11 ±0.04 

-10.43 ±0.21 

39.89 

(2) 

NGC3049 

9.71 ±0.19 

-0.41 ± 0.06 

-10.12 ±0.20 

40.58 

(2) 

NGC3077 

9.46 ± 0.20 

-1.06 ±0.02 

-10.52 ±0.20 

39.98 

(4) 

M8IdwB 

7.24 ± 0.32 

-2.36 ±0.34 

-9.61 ±0.47 

38.37“ 

(4) 

NGC3I90 

10.74 ±0.20 

-0.47 ±0.18 

-11.21 ±0.27 

39.59 

(2) 

NGC3I84 

10.49 ±0.21 

-0.07 ± 0.09 

-10.56 ±0.23 

40.93 

(2) 

NGC3198 

10.35 ±0.21 

-0.12 ±0.10 

-10.47 ±0.23 

40.87 

(2) 

IC2574 

8.99 ± 0.20 

-1.18 ±0.08 

-10.17 ±0.22 

40.01 

(2) 

NGC3265 

9.55 ±0.21 

-0.65 ± 0.09 

-10.21 ±0.23 

40.22 

(2) 

NGC335I 

10.47 ±0.20 

-0.27 ± 0.08 

-10.74 ±0.22 

40.58 

(2) 

NGC352I 

10.99 ± 0.20 

0.37 ± 0.05 

-10.62 ±0.21 

41.17 

(2) 

NGC362I 

10.20 ±0.20 

0.01 ± 0.05 

-10.19 ±0.20 

41.11 

(2) 

NGC3627 

10.79 ±0.20 

0.29 ± 0.04 

-10.50 ±0.20 

41.11 

(2) 

NGC3773 

9.14 ±0.20 

-0.93 ± 0.02 

-10.08 ±0.20 

40.21 

(2) 

NGC3938 

10.60 ±0.19 

0.21 ± 0.02 

-10.39 ±0.19 

41.20 

(2) 

NGC4236 

9.08 ±0.19 

-0.87 ± 0.06 

-9.95 ± 0.20 

40.33 

(4) 

NGC4254 

10.73 ±0.20 

0.49 ± 0.04 

-10.25 ±0.21 

41.38 

(2) 

NGC432I 

10.86 ±0.20 

0.35 ± 0.04 

-10.51 ±0.21 

41.20 

(2) 

NGC4536 

10.44 ±0.19 

0.17 ±0.11 

-10.26 ±0.22 

40.88 

(2) 

NGC4559 

9.82 ±0.21 

-0.40 ± 0.05 

-10.23 ±0.21 

40.71 

(2) 

NGC4569 

10.51 ±0.19 

-0.45 ± 0.08 

-10.96 ±0.21 

40.29 

(2) 

NGC4579 

11.03 ±0.20 

-0.05 ± 0.06 

-11.09 ±0.21 

40.86 

(2) 

NGC4594 

11.19 ±0.20 

-0.75 ±0.17 

-11.95 ±0.26 

39.70 

(4) 

NGC4625 

9.27 ±0.19 

-1.23 ±0.15 

-10.51 ±0.24 

39.81 

(2) 

NGC463I 

10.44 ±0.20 

0.27 ± 0.04 

-10.17 ±0.20 

41.19 

(2) 

NGC4725 

10.85 ±0.20 

-0.09 ± 0.20 

-10.94 ±0.28 

40.98 

(13) 

NGC4736 

10.51 ±0.20 

-0.36 ±0.10 

-10.87 ±0.22 

40.48 

(2) 

DDOI54 

7.35 ±0.41 

-2.24 ± 0.20 

-9.59 ± 0.46 

38.57“ 

(2) 

NGC4826 

10.48 ±0.20 

-0.57 ±0.10 

-11.05 ±0.22 

40.18 

(2) 

DDOI65 

8.01 ±0.20 

-2.30 ± 0.26 

-10.31 ±0.33 

38.48“ 

(2) 

NGC5055 

10.77 ±0.20 

0.11 ±0.06 

-10.66 ±0.21 

40.91 

(2) 

NGC5398 

8.99 ± 0.02 

-1.18 ±0.06 

-10.18 ±0.07 

39.94 

(13) 

NGC5408 

8.60 ± 0.20 

-1.06 ±0.02 

-9.67 ± 0.20 

40.15 

(2) 

NGC5457 

10.67 ±0.20 

0.32 ±0.12 

-10.35 ±0.23 

41.33 

(4) 

NGC5474 

9.17 ±0.17 

-1.11 ±0.25 

-10.28 ±0.30 

40.01 

(2) 

NGC5713 

10.50 ±0.22 

0.32 ±0.10 

-10.18 ±0.24 

40.96 

(2) 

NGC5866 

10.80 ±0.20 

-0.38 ±0.13 

-1L18±0.24 

40.53 

(14) 

NGC6946 

10.77 ± 0.20 

0.53 ± 0.04 

-10.24 ±0.21 

41.53 

(2) 

NGC733I 

11.15 ±0.20 

0.49 ± 0.07 

-10.66 ±0.21 

41.22 

(2) 

NGC7793 

9.72 ±0.19 

-0.57 ± 0.05 

-10.28 ±0.20 

40.56 

(2) 


Notes: (1) Log of stellar m asses i n LogiMo], estimated with the formula from |Eskew et al.|f2012 1; (2) Log of star-formation rates from Hct+Ltir, in 
LogiMg yr“'] (see Section 2.6.1 1 ; (3) Log of specific star-formation rates defined by SFR/Ms,aj, in Log[yr“']; (4) Log of Hn luminosities corrected 
for underlying stellar absorption, Nil line contamination and foreground Galactic extinction, in Log[erg/s]; (5) References for La„. 

“ SFR estimated from |Lee et al.| ( [^09| l (see Section [2.6.11 . 


References for Ha luminosities : 

(l)|Schmitt et al.|l|2006 

; (2)|Kennicutt et al. (|20091l ; (3) 

Moustakas & Kennicutt|(|2006J ; (4)|Kennicutt et al. 

(|2008) ; (5)|Gil de Paz et al.|j2003 

; (6)|Popescu & Hopp 

12000J ; (7)|Iglesias-Paramo 8i Vilchez|(|1997|l ; (8)|Pustilnik et al.| 

2004) ; (9)|Izotov 

let al.|(|1994|l ; (10)|Izotov & Thuanj 

1998|l ; (1 l)|Terlevich et al.|(119911l ; (12)|Moles et al.|(|1994|l ; (13)|Skibba et al.|(|201 IJ ; (14) 

Plana etal.|11998|l 
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Table 5. Updated Herschel flux densities for the DGS sample. 


Name 

F 70 

[Jy] 

PACS 

Fioo 

[Jy] 

F 160 

[Jy] 

Map 

method" 

F 250 

[Jy] 

SPIRE 

F 350 

[Jy] 

F 500 

[Jy] 

Haro 11 

6.34 ± 0.32 

5.04 ± 0.25 

2.15 ± 0.11 

S 

0.59 ± 0.05 

0.22 ± 0.02 

0.080 ± 0.009 

Haro2 

5.09 ± 0.25 

5.47 ± 0.27 

3.66 ±0.18 

S 

1.23 ± 0.10 

0.52 ± 0.04 

0.17 ±0.02 

Haro3 

5.59 ± 0.28 

6.09 ± 0.30 

4.73 ± 0.24 

S 

1.72 ±0.14 

0.76 ± 0.06 

0.25 ± 0.02 

He2-10 

26.9 ± 1.3 

26.4 ± 1.3 

17.7 ± 0.9 

S 

6.47 ± 0.52 

2.56 ±0.21 

0.90 ± 0.08 

HS0017+1055 

0.046 ± 0.005 

0.033 ± 0.004 

0.019 ± 0.004 

P 

< 0.030 

< 0.050 

< 0.045 

HS0052-(-2536 

0.22 ± 0.01 

0.21 ± 0.01 

0.139 ±0.008 

P 

0.056 ± 0.007 

0.030 ± 0.010 

0.021 ± 0.008 

HS0822-I-3542 

< 0.041* 

< 0.048* 

0.034 ± 0.003 

P 


- 

- 

HS1222-I-3741 

0.025 ± 0.004 

< 0.036 

< 0.022 

P 


- 

- 

HS1236+3937 

< 0.029 

< 0.035 

< 0.028 

P 

< 0.030 

< 0.050 

< 0.045 

HS1304-E3529 

0.121 ±0.007 

0.150 ±0.009 

0.069 ± 0.005 

P 

0.029 ± 0.005 

< 0.050 

< 0.045 

HS1319-H3224 

0.012 ± 0.003 

0.013 ± 0.002 

< 0.015 

P 


- 

- 

HS1330+3651 

0.093 ± 0.006 

0.112 ±0.007 

0.091 ± 0.005 

P 


- 

- 

HS1442-h4250 

0.09 ± 0.01 

< 0.054 

< 0.047 

P 

< 0.030 

< 0.050 

< 0.045 

HS2352-H2733 

0.039 ± 0.003 

0.016 ± 0.002 

< 0.016 

P 

< 0.030 

< 0.050 

< 0.045 

IZwl8 

0.045 ± 0.003 

0.018 ± 0.002 

<0.011 

P 

< 0.030 

< 0.050 

< 0.045 

ICIO 

145. ±7. 

211. ± 11. 

209. ± 10. 

s 

99.5 ± 30.9 

48.1 ± 14.9 

19.3 ± 6.0 

IIZw40 

6.66 ± 0.33 

6.01 ± 0.30 

3.27 ±0.16 

s 

1.27 ± 0.11 

0.53 ± 0.05 

0.15 ±0.01 

Mrkl089 

4.86 ± 0.24 

5.17 ± 0.26 

4.41 ± 0.22 

s 

1.72 ±0.14 

0.79 ± 0.07 

0.29 ± 0.03 

Mrkl450 

0.30 ± 0.02 

0.25 ± 0.01 

0.127 ±0.007 

p 

0.046 ± 0.006 

< 0.050 

< 0.045 

Mrkl53 

0.28 ± 0.01 

0.30 ± 0.02 

0.137 ±0.009 

p 

0.045 ± 0.007 

< 0.050 

< 0.045 

Mrk209 

0.36 ± 0.02 

0.37 ± 0.02 

0.17 ±0.01 

s 

0.059 ± 0.007 

0.033 ± 0.009 

< 0.045 

Mrk930 

1.20 ±0.06 

1.48 ± 0.07 

0.96 ± 0.05 

s 

0.39 ± 0.03 

0.19 ±0.02 

0.08 ± 0.01 

NGC1140 

4.24 ± 0.21 

4.44 ± 0.22 

4.47 ± 0.22 

s 

1.89 ±0.16 

0.92 ± 0.08 

0.33 ± 0.04 

NGC1569 

61.7 ±3.1 

58.9 ± 2.9 

37.4 ± 1.9 

s 

12.1 ± 1.0 

5.03 ±0.41 

1.85 ±0.16 

NGC1705 

1.10 ±0.06 

1.37 ± 0.07 

1.14 ±0.06 

S" 

0.59 ± 0.05 

0.30 ± 0.03 

0.12 ±0.02 

NGC2366 

5.82 ± 0.29 

7.05 ± 0.35 

4.30 ± 0.22 

s 

2.00 ± 0.17 

1.07 ±0.09 

0.45 ± 0.04 

NGC4214 

26.8 ± 1.3 

35.0 ± 1.8 

33.0 ± 1.6 

s 

18.6 ± 1.5 

9.96 ±0.81 

4.53 ± 0.37 

NGC4449 

50.4 ± 2.5 

79.6 ± 4.0 

75.6 ±3.8 

s 

32.2 ± 2.6 

14.8 ± 1.2 

5.90 ± 0.48 

NGC4861 

2.31 ±0.12 

2.27 ±0.11 

2.36 ±0.12 

s 

1.06 ±0.09 

0.55 ± 0.05 

0.23 ± 0.03 

NGC5253 

33.7 ± 1.7 

33.5 ± 1.7 

22.5 ± 1.1 

s 

7.78 ± 0.63 

3.59 ± 0.29 

1.37 ±0.11 

NGC625 

6.14 ±0.31 

9.98 ± 0.50 

8.32 ± 0.42 

s 

4.34 ± 0.35 

2.18 ±0.18 

0.95 ± 0.08 

NGC6822 

56.2 ± 2.9 

72.5 ±3.7 

79.0 ± 4.0 

s 

49.7 ± 15.5 

30.8 ± 9.6 

16.1 ±5.0 

Pox186 

0.051 ±0.003 

0.059 ± 0.005 

0.067 ± 0.001 

p 

0.042 ± 0.006 

< 0.050 

< 0.045 

SBS0335-052 

0.056 ± 0.004 

0.024 ± 0.001 

0.007 ± 0.001 

p 

< 0.030 

< 0.050 

< 0.045 

SBS1159-P545 

0.019 ± 0.003 

0.019 ± 0.003 

<0.018 

p 

< 0.030 

< 0.050 

< 0.045 

SBS1211-P540 

0.034 ± 0.003 

0.018 ± 0.002 

0.013 ± 0.002 

p 

< 0.030 

< 0.050 

< 0.045 

SBS1249-P493 

0.032 ± 0.005 

< 0.034 

< 0.042 

p 

< 0.030 

< 0.050 

< 0.045 

SBS1415-P437 

0.18 ±0.01 

0.16 ±0.01 

0.065 ± 0.007 

p 


- 

- 

SBS1533-P574 

0.23 ± 0.01 

0.30 ± 0.02 

0.174 ±0.010 

s 

<0.145 

< 0.050 

< 0.045 

Tol0618-402 

< 0.014 

<0.011 

<0.013 

p 


- 

- 

Tol 1214-277 

0.017 ± 0.003 

0.018 ± 0.002 

<0.018 

p 

< 0.030 

< 0.050 

< 0.045 

UGC4483 

0.11 ±0.02 

0.27 ± 0.02 

0.091 ±0.010 

S" 

0.015 ± 0.005 

< 0.050 

< 0.045 

UGCA20 

< 0.052 

< 0.057 

< 0.048 

p 


- 

- 

UM133 

0.12 ±0.02 

0.13 ±0.02 

0.078 ± 0.009 

prf 

0.025 ± 0.006 

< 0.050 

< 0.045 

UM311 

3.13 ±0.16 

5.52 ± 0.28 

6.12 ±0.31 

s 

3.79 ±0.31 

1.91 ±0.16 

0.81 ±0.07 

UM448 

5.23 ± 0.26 

- 

3.32 ±0.17 

s 

0.95 ± 0.08 

0.36 ± 0.03 

0.11 ±0.01 

UM461 

0.21 ±0.01 

0.145 ± 0.009 

0.113 ±0.007 

p 

0.025 ± 0.006 

0.024 ± 0.010 

< 0.045 

VIIZw403 

0.47 ± 0.03 

0.56 ± 0.03 

0.34 ± 0.02 

p 

0.14 ±0.01 

0.062 ± 0.008 

0.027 ± 0.008 


Notes : 

“ P = PhOtPrOJECT, S = SCANAMORPHOS. 

* These upper limits come from the SED modelling (see Appendix]^. 

For these galaxies, the PACS map-making method was changed from PhotProject to Scanamorphos since 
For these galaxies, the PACS map-making method was changed from Scanamorphos to PhotProject since 


Remy^^yeret^ 1 2013 . 

Remy-Ruyer et al. 2013 i. 
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Table 6. Spitzer IRAC flux densities for the DGS sample. 


Source 

F 3.6 

[mJy] 

F 4.5 

[mJy] 

F 5.8 

[mJy] 

Fg.o 

[mJy] 

Aperture radius 
["] 

Aper corr ? 

Haro 11 

22.4 ± 3.2 

32.3 ± 4.6 

77.1 ± 10.9 

164. ± 23. 

45“ 

yes 

Haro2 

24.7 ± 3.5 

18.2 ± 2.6 

49.3 ± 7.0 

257. + 36. 

50“ 

yes 

Haro3 

28.2 ± 4.0 

21.5 + 3.1 

43.9 ± 6.3 

116. ± 16. 

60“ 

yes 

He2-10 

134. ± 19. 

99.2+ 14.1 

296. ± 42. 

836.+ 118. 

108“ 

yes 

HS0017+1055 

0.17 ± 0.02 

0.31 + 0.03 

0.76 + 0.08 

1.53 + 0.16 

12* 


HS0052+2536 

0.7 ±0.1 

0.32 + 0.09 

0.7 ± 0.2 

3.17 + 0.46 

17“ 

yes 

HS0822+3542 

0.12 + 0.01 

0.10 + 0.01 

0.10 + 0.02 

0.09 + 0.02 

5* 

yes 

HS1222+3741 

0.25 ± 0.03 

0.19 + 0.02 

0.27 ± 0.05 

0.57 ± 0.07 

14“ 


HS1236+3937 

<0.56 

<0.38 

<0.66 

<0.32 

15“ 


HS1304+3529 

0.38 ± 0.08 

0.41 ± 0.09 

<0.78 

0.6 ±0.1 

18“ 


HS 1319+3224 

0.061 ± 0.010 

0.051 + 0.008 

0.10 + 0.02 

0.12 + 0.02 

8“ 


HS 1330+3651 

0.56 ± 0.09 

0.5+ 0.1 

<0.80 

1.0 + 0.2 

20“ 


HS 1442+4250 

- 

1.24 + 0.21 

- 

1.59 ± 0.25 

51“ 

yes 

HS2352+2733 

<0.04 

<0.04 

<0.13 

<0.11 

15“ 


IZwlS 

0.36 ± 0.05 

0.35 + 0.05 

0.39 ± 0.07 

0.53 ± 0.08 

12* 

yes 

ICIO^ 

- 

- 

- 

- 

- 

- 

IIZw40 

18.1+4.3 

19.9 ± 3.5 

39.8 + 5.7 

105. + 15. 

45 X 33“ 

yes 

Mrkl089 

20.4 ± 2.9 

15.5 + 2.2 

31.4 + 4.4 

82.2+ 11.6 

75“ 

yes 

Mrkl450 

1.22 + 0.17 

1.01 ± 0.14 

0.9 ±0.1 

2.06 ± 0.29 

20“ 

yes 

Mrkl53 

2.29 ± 0.33 

1.70 + 0.24 

1.76 + 0.25 

2.75 + 0.39 

25* 

yes 

Mrk209 

2.90 ± 0.41 

2.21 ± 0.32 

1.78 + 0.27 

2.20 ± 0.32 

39“' 

yes 

Mrk930 

2.66 ± 0.38 

2.73 + 0.39 

4.17 + 0.59 

10. + 1. 

13* 

yes 

NGC1140 

32.8 + 4.8 

23.3 ± 3.4 

41.0 + 5.8 

93.9+ 13.3 

118“ 

yes 

NGC1569 

312. ±45. 

244. ± 35. 

324. ± 46. 

533. ±75. 

150“ 

yes 

NGC1705 

25.5 ± 5.9 

18.3 ± 4.2 

15.0 + 2.8 

17.5 + 2.6 

72“ 

yes 

NGC2366 

61.5 ± 12.4 

45.7 ± 10.4 

35.8 + 7.6 

52.8 + 8.1 

200“^ 

yes 

NGC4214 

334. ± 47. 

253. ± 36. 

280. ± 40. 

748. ± 106. 

205 X 160“ 

yes 

NGC4449 

472. ± 67. 

342. + 48. 

709. + 100. 

1622. + 229. 

190 X 170“ 

yes 

NGC4861 

18.3 + 2.6 

13.8 + 2.0 

16.8 + 2.5 

17.7 ± 2.5 

100 X 50“ 

yes 

NGC5253 

235. ± 34. 

252. ± 36. 

458. ± 65. 

812. ± 115. 

120“ 

yes 

NGC625 

110. ± 16. 

79.7 ±11.4 

71.3 ± 10.1 

147. ±21. 

170“ 

yes 

NGC6822 

1780. ± 272. 

1240. ± 196. 

951. ± 141. 

1226. + 175. 

440“ 

yes 

Pox186 

0.18 + 0.02 

0.18 + 0.02 

0.14 + 0.03 

0.64 ± 0.07 

9* 


SBS0335-052 

0.78 + 0.08 

1.64 + 0.17 

4.50 ± 0.45 

12.3 + 1.2 

10“ 


SBS1159+545 

0.11+0.01 

0.10 + 0.01 

0.12 + 0.03 

0.59 ± 0.07 

12“' 


SBS1211+540 

0.16 + 0.02 

0.08 ± 0.01 

<0.12 

0.30 ± 0.05 

15“ 


SBS 1249+493 

0.08 ± 0.02 

0.07 + 0.02 

<0.20 

0.26 ± 0.04 

5* 


SBS1415+437 

- 

0.8 ±0.1 

- 

0.50 ± 0.08 

34“ 

yes 

SBS1533+574 

2.09 ± 0.34 

1.52 ± 0.24 

2.26 ± 0.33 

3.32 ± 0.48 

30“ 

yes 

Tol0618-402 

2.68 + 0.38 

1.64 + 0.23 

1.13 + 0.16 

0.9+ 0.1 

18“ 

yes 

ToH214-277 

0.076 ± 0.009 

0.09 ± 0.01 

0.15 + 0.03 

0.25 ± 0.03 

8* 


UGC4483 

1.69 + 0.25 

1.02 + 0.15 

<0.62 

0.9 ±0.1 

43 X 24“ 

yes 

UGCA20 

<0.85 

<0.60 

< 1.07 

<0.46 

20“ 


UM133 

0.39 ± 0.08 

0.28 ± 0.06 

<0.46 

<0.27 

8“ 


UM311 

- 

33.4 ± 4.7 

- 

114. ± 16. 

115“ 

yes 

UM448 

17.7 ± 2.5 

13.2+ 1.9 

34.1+4.8 

83.0+ 11.7 

53x31“ 

yes 

UM461 

0.8 ±0.1 

0.61 ± 0.09 

0.8 ±0.1 

1.58 + 0.29 

17“ 

yes 

VIIZw403 

2.86 + 0.41 

2.23 ± 0.32 

1.69 + 0.25 

2.56 ± 0.37 

40“ 

yes 


Notes: 

“ The aperture is the same as the one used for Herschel. 

* The Herschel aperture has been shorten to avoid a contaminating source. 

^ The Herschel aperture has been adapted to match the peculiar morphology of the source in the NIR. 

The Herschel aperture has been enlarged to encompass all of the NIR emission. 

'■ The IRAC map is smaller than the Herschel aperture. The aperture thus had to be shorten, 
f The IRAC maps do not cover the full galaxy. Thus we do not report flux densities for this source. 
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Table 7. WISE flux densities for the DGS sample. 


Source 

F3.4 

[mjy] 

F4.6 

[mJy] 

Fi2 

[mJy] 

F22 

[mJy] 

Aperture radius 
["] 

Haro 11 

17.9 ± 1.9 

32.6 ±3.4 

353. ± 19. 

2054. ± 148. 

45“, 90" 

Haro2 

24.5 ± 2.6 

16.5 ± 1.7 

152. ± 8 . 

760. ± 55. 

50“, 65*’ 

Haro3 

25.8 ± 2.7 

17.7 ± 1.9 

162. ± 9. 

766. ± 56. 

60“ 

He2-10 

133.± 17. 

92.7 ± 11.2 

996. ± 53. 

5018. ± 362. 

81* 

HS0017-H1055 

0.16 ±0.03 

0.25 ± 0.04 

3.83 ± 0.49 

11.7 ± 2.0 

12“, 25* 

HS0052-H2536 

0.65 ± 0.08 

0.52 ± 0.07 

3.71 ± 0.40 

16.2 ±2.3 

17“, 23“ 

HS0822-H3542 

0.106 ±0.009" 

0.082 ± 0.008" 

< 0.47" 

< 3.85" 

- 

HS1222-H3741 

0.15 ±0.03 

0.15 ±0.03 

1.0 ± 0 . 2 " 

5.94 ± 0.73" 

14“, - 

HS1236-H3937 

0.22 ± 0.05 

0.16 ± 0.04 

<0.65 

< 2.72 

15“ 

HS1304-H3529 

0.4 ±0.1 

0.43 ± 0.09 

1.78 ±0.21" 

10.5 ± 0.9" 

18“, - 

HS1319-H3224 

0.06 ± 0.02 

0.07 ± 0.02 

< 0.70" 

< 1.93" 

12 *, - 

HS1330-H3651 

0.51 ±0.07 

0.36 ± 0.07 

1.23 ± 0.26 

5.31 ± 1.26 

20 “ 

HS1442-H4250 

1.20 ±0.14 

0.8 ± 0.1 

1.54 ± 0.33 

<9.40 

32* 

HS2352-H2733 

0.15 ±0.04 

0.13 ±0.04 

< 1.03" 

< 2.98" 

15“,- 

IZwl 8 

0.33 ± 0.07 

0.39 ± 0.07 

0.8 ± 0.3 

4.64 ± 1.19 

23* 

ICIO 

1454. ± 154. 

917. ±98. 

3521. ± 191. 

9991. ± 728. 

205 X 150*, 230* 

IIZw40 

18.6 ± 2.0 

19.3 ±2.1 

335. ± 18. 

1540. ± 112. 

45 X 33“, 90* 

Mrkl089 

19.0 ± 2.3 

13.9 ±3.1 

102 . ± 6 . 

453. ± 34. 

75“ 

Mrkl450 

1.10 ± 0.12 

0.73 ± 0.09 

8.00 ± 0.74 

49.3 ± 4.8 

20“, 40* 

Mrkl53 

2.03 ± 0.23 

1.34 ± 0.17 

4.61 ± 0.41 

28.1 ±3.9 

25“ 

Mrk209 

2.61 ±0.39 

1.70 ± 0.31 

7.30 ±0.85 

46.4 ± 5.4 

39“ 

Mrk930 

2.61 ±0.30 

2.26 ± 0.27 

25.8 ± 2.6 

170. ± 14. 

34*, 60“ 

NGC1140 

32.7 ± 3.5 

21.6 ±2.3 

94.9 ± 5.5 

367. ± 28. 

118“ 

NGC1569 

359. ± 24. 

240. ± 47. 

1021 . ± 116. 

6971. ± 503. 

120 * 

NGC1705 

25.7 ± 2.7 

16.1 ± 1.7 

20.3 ± 1.5 

46.2 ±6.1 

75 X 46*, 72“ 

NGC2366 

75.8 ±5.7 

39.4 ± 6.2 

121 . ± 12 . 

612. ±65. 

260 X 125*, 150“ 

NGC4214 

294. ± 31. 

184. ± 19. 

546. ± 29. 

1830. ± 133. 

181“ 

NGC4449 

483. ± 50. 

307. ± 32. 

1323.8 ±70.0 

3119. ± 227. 

190 X 170“, 250“ 

NGC4861 

18.1 ± 1.1 

10.1 ± 1.9 

52.8 ± 5.7 

330. ± 26. 

120 X 50*, 120“ 

NGC5253 

212 . ± 22 . 

255. ± 27. 

2002. ± 105. 

10345. ± 746. 

120 X 80*, 120“ 

NGC625 

106. ± 11 . 

64.7 ± 6.8 

194. ± 11. 

797. ± 59. 

200 X 75 

NGC6822‘' 

- 

- 

- 

- 

- 

Pox186 

<0.34 

<0.44 

2.67 ± 0.33 

11.2 ± 1.4 

25* 

SBS0335-052 

0.45 ± 0.05 

1.45 ± 0.16 

23.4 ± 1.4 

67.8 ±5.8 

14*, 35* 

SBSl 159-1-545 

<0.09 

< 0.11 

1.0 ±0.3 

4.44 ± 1.38 

12“, 25* 

SBS1211-t540 

0.133 ±0.009 

0.10 ±0.03 

<0.80 

<3.60 

15“ 

SBS12494-493 

< 0.05" 

< 0.07" 

< 1.17 

<3.76 

-, 20 * 

SBS1415-h437 

0.8 ± 0.1 

0.53 ± 0.08 

2.30 ± 0.55 

15.5 ±2.8 

28 X 15*, 34“ 

SBS1533-h574 

1.52 ±0.19 

1.28 ± 0.17 

9.76 ± 0.84 

53.3 ± 4.7 

30“, 50* 

T 0 IO 61 8-402 

2.33 ± 0.33 

1.32 ±0.21 

<0.71 

< 1.73 

18“ 

Tol 1214-277 

<0.13 

<0.14 

0.8 ± 0 . 2 " 

< 6.89" 

12 “, - 

UGC4483 

1.12 ±0.26 

< 1.01 

<2.30 

< 12.72 

43 X 24 

UGCA20 

0.26 ± 0.08 

<0.24 

<0.74 

< 3.27 

20 “ 

UM133 

< 10.83 

< 6.74 

< 1.49 

<6.70 

26“ 

UM311 

48.5 ± 3.0 

29.1 ±3.5 

120. ± 7. 

292. ± 25. 

115“ 

UM448 

15.4 ± 1.6 

11.7 ± 1.3 

103. ± 6 . 

552. ± 40. 

64“ 

UM461 

0.58 ± 0.07 

0.52 ± 0.07 

5.98 ± 1.02 

29.4 ± 4.9 

17“, 35* 

VIIZw403 

2.59 ± 0.29 

1.89 ± 0.23 

4.54 ± 0.61 

26.3 ± 3.6 

40“ 


Notes: 

In the aperture column, a single value indicates that only one aperture has been used for the four wavelengths. Two values separated by a comma 
indicate that a different aperture has been used for WISEl and WISE2 on one side (first value), and WISE3 and WISE4 on the other side 
(second value, see text for details). For several galaxies elliptical apertures were used to match the morphology of the source and are indicated by 
semi-major axis x semi-minor axis. A dash indicates that we take the profile fit photometry from the All WISE database. 

“ The aperture is the same as that used for IRAC or Herschel photometry. 

* The aperture has been adapted to match the peculiar morphology of the source in the NIR/MIR or to avoid a contaminating source. 

^ This galaxy is not resolved at these wavelengths and we thus use the profile fit photometry provided by the WISE database. 

The WISE map does not cover the full galaxy. Thus we do not report flux densities for this source. 



Table 8. 2MASS and IRAS flux densities from the literature for the DGS sample. 
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Table 9. (continued) DGS and KINGFISH dust parameters. 
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Appendix A: Update of the Herschel data 

A.1. PACS 


In [Remy-Ruyer et aL (2013 i, the PACS data for the DGS galax¬ 
ies were processed either with PhotProject (provided by HIPE) 
for point sources or with Scanamorphos (version 14, Roussel 
2013|l for extended sources. The update of Scanamorphos to 


try on the new data, depending if the source is extended or point 
like. For the extended sources, there are two steps to apply before 
performing aperture photometry: convert the data from point- 
source calibration to extended-source calibration, via the Kp,oe 
factof^ and convert the data from Jy/beam to Jy/pix, using the 
beam area. The values for the KptoE factor, given the SPIRE Data 
Reduction Guid^ are 0.99858, 1.00151 and 0.99933 at 250, 
350 and 500 pm respectively. For the “slightly extended” sources 


This factor was denoted 
The SPIRE Data Reduction 


m 


Remy-Ruyer et al. 1 2013[l. 


version 23 implies a significant update on the red band at 160 
pm due to the PACS distortion flat-held that is now taken into 
account. This effect was already included in PhotProject. For 
the DGS data, the Scanamorphos maps are thus reprocessed 
with v23 of the algorithm for the three PACS bands for con¬ 
sistency. The photometry is done with the same method as in 
Remy-Ruyer et al. ( 2013| ), and using the exact same apertures. 
There are three changes to note : 

- NGC1705 - The vl4 maps of Scanamorphos for this galaxy 
presented a non-uniform background which compelled us to 
choose the PhotProject maps in Remy-Ruyer et al. ( |2013| ) (see 
note 1 in their Table 2). This has been corrected in the new ver¬ 
sion of the maps and so we now consider the Scanamorphos 
maps for this galaxy. 

- UGC 4483 - The new Scanamorphos maps now manage to 
reconstruct better the emission for this galaxy at 160 pm yielding 
a strong detection (9cr) in the new map. Moreover, this detection 
at 160 pm is more coherent with the rest of the FIR/submm pho¬ 
tometry (detection at 250 pm by SPIRE). So we now consider 
the Scanamorphos maps for this galaxy and not the PhotProject 
maps as in Remy-Ruyer et al. (2013 1 . 

- UM 133 - For this galaxy Remy-Ruyer et al. ( 2013| l noted 
that the 70 pm flux density might present some discrepancies 
with others FIR measurements. Using PhotProject maps 
instead of Scanamorphos maps yields a 70 pm flux density 
which is more consistent with the FIR shape of the SED. So we 
now consider the PhotProject maps for this galaxy. 

The new flux densities for the DGS are presented in Table 
1^ For the KINGFISH PACS data, we applied a corrective factor 
of 0.925 to the 160 pm flux densities of [Dale et al.| ( |2012] l (H. 
Roussel, priv. com.). 


A.2. SPIRE 

There have been two different updates for the SPIRE data : the 
calibration and the beam model used to compute the beam ar- 

ea^3 

The SPIRE maps for the DGS have been reprocessed with 
HIPE vl2, as in Remy-Ruyer et al. ([2013 1 . We apply the same 
method as in Remy-Ruyer et al. (2013 1 to redo the photome¬ 


(see Remy-Ruyer et al.|2013 i we do not apply this KptoE factor. 
The beam areas also depend on the spectral shape of the consid¬ 
ered source, which means that the beam areas will be different 
for each galaxy. To produce the final maps, the SPIRE pipeline 
assumes a spectral shape with a dependence Fy oc v“*, which 
corresponds to beam areas of 465, 823 and 1769 arcsec^ at 250 
pm, 350 pm and 500 pm respectively. These are the beam areas 
we used for the second step (identical to those used in [Remy- 
Ruyer et al.|2013| l, before performing aperture photometry. 

To account for the different FIR spectral shapes and thus 
beam areas, we need to apply a colour correction factor. The 
colour correction factors, Kcoi, have been tabulated in the SPIRE 
Data Reduction Guide for spectral shapes with Fy cc •/’. For 
point sources and extended sources, the Kcoi are given in Table 
6.13 of the SPIRE Data Reduction Guide. For slightly extended 
sources, we apply the point source colour correction times the 
effective beam area ratios tabulated in Table 6.8 of the SPIRE 
Data Reduction Guide. Note that this effective beam area ratio is 
already taken into account in the colour correction for extended 
sources. The method we use to colour correct our flux densities 
is the following: we fit a line in logarithmic space to the three 
SPIRE flux densities to find the spectral slope a. We include 
PACS 160 pm data if there are some SPIRE non-detections. 
Once a is determined, we apply the corresponding colour cor¬ 
rection factors to the flux densities. We iterate these few steps 
until a does not vary by more than 0.1%. Usually two iterations 
are enough. 

In Table we give the updated fluxes densities before the 
colour correction step. On the other hand, the colour-corrected 
flux densities are the ones used in the fitting procedure. The 
colour corrections are of the order of 6%, 4% and 4% at 250 
pm, 350 pm, and 500 pm respectively for the DGS. 


The beam areas are needed to convert the original data in Jy/beam 
to Jy/pix for the aperture photometry. 


Guide, v3.0, is available at: 
http://herschel.esac.esa.int/hcss-doc-l 2.0/print/spire _drg/spire_drg.pdf 
#spire_drg. 


The SPIRE maps of [Dale ef^ ( |20T2| l for the KINGFISH 
sample were reduced with HIPE v5, and aperture photometry 
was performed on the maps to extract the flux densities. Since 
then, there have been two updates of the calibration with the re¬ 
leases of v7 and vl 1 of HIPE. To account for the update at v7 we 
multiply the flux densities from Dale et al. ( 2012[ ) by 1.0, 1.0067 
and 1 .C p^ at 250, 350 and 500 pm respectively. To account for 
vl 1 update, we multiply the flux densities by 1.0253, 1.0250 and 
1.0125^^ at 250, 350 and 500 pm respectively. No KptoE factor 
had previously been applied to the flux densities (D. Dale, priv. 
com.). Assuming that the KINGFISH sources are all extended, 
we also apply the Kp,oe factor to all of the galaxies. We also 
convert the fluxes to match the beam areas mentioned previously 
(465, 823 and 1769 arcsec^). Then these updated flux densities 
are colour corrected following the method described for the DGS 
galaxies above. The colour corrections are of the order of ~ 5% 
for the three SPIRE wavelengths for the KINGFISH galaxies. 


Appendix B: DGS IRAC and IRS Observing Logs 

Appendix C: Quality checks for IRAC and WISE 
DGS photometry 

C.1. IRAC: Comparison to previous literature measurements 

We compare IRAC flux densities previously available in the liter¬ 
ature and the flux densities derived in this work. We use the ratios 
of our IRAC measurement to the literature IRAC flux densities; 
a ratio of 1 cotTesponds to a very good agreement. The compari- 


G. J. Bendo, priv. com. 
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Fig. C.l. (left panels) Comparison of our IRAC flux densities and literature IRAC flux densities: F(Ao)IFLIT(Ao) flux density 
ratios as a function of our IRAC flux density, F{Aq), at 3.6 pm, 4.5 pm, 5.8 pm, and 8.0 pm. The average uncertainties are 17%, 
17%, 19% and 17% at 3.6, 4.5, 5.8 and 8.0 pm. Colours distinguish the reference for the literature measurement, (right panels) 
Ratios of WISE flux densities to other MIR measurements (either IRAC, IRAS or MIPS) as a function of WISE flux density at 3.4 
pm, 4.6 pm, 12 pm and 22 pm. The average uncertainties are 20%, 22%, 15% and 12% at 3.4, 4.6, 12 and 22 pm respectively. On 
all panels, we show the unit ratio as a solid line as well as the average uncertainties on the ratio in each bands as dashed lines. 


son is presented in Eig. |C.l| and shows good agreement between 
the two measurements, even if some outliers are present: 


- Haro 2 - The comparison to the literature values from Dale 


et al. (2007 1 agrees for the first three IRAC bands (ratios of 0.92, 


0.96 and 0.93 at 3.6 pm, 4.5 pm, and 5.8 pm respectively) but 
not for the IRAC 8.0 pm band where the ratio goes to ~ 2. 


Recomputing the flux density with the aperture given in Dale 
|et al.| ( |2007| l gives a flux in very close agreement to the one 
we find with our own aperture (ratio of 1.009 at 8.0 pm), and 
also for the three other IRAC bands. Additionally, we find that 
[Marble et al.| ( |2010j l performed photometry for Haro 2 in an IRS- 
matched aperture and their result is in good agreement with ours 
(ratio = 1.036 at 8.0 pm). 

- NGC 4861 - UGC 4483 - Tol 1214-277 - Our measurements 


are systematically larger than the flux densities from 

Engelbracht 

et al. ([2008 1. No precise information is given in 

Engelbracht 

et al. ( 2008|l about the size of the apertures used for IRAC pho- 


tometry, so it is difficult to asses that we are really compar¬ 
ing similar measurements. Nonetheless, a possible explanation 
for NGC 4861 could be that we include Mrk 39 with which 
NGC 4861 is interacting in our aperture and that [Engelbrachtj 
jet al.| ( |2008l l did not. 

Excluding points where one of the two measurements (ours 
or literature values) is an upper limit, we get a mean ratio of our 
IRAC measurements to those in the literature of 1.05 +0.11 for 
IRAC 3.6 pm, 1.06 + 0.11 for IRAC 4.5 pm, 1.02 + 0.14 for 
IRAC 5.8 pm and 1.01 + 0.11 for IRAC 8.0 pm. This is to be 
compared to an average uncertainty on the ratios of 17%, 17%, 
19% and 17% at 3.6, 4.5, 5.8 and 8.0 pm. The error on the mean 
ratio is lower than the average uncertainty on the ratios for all of 
the bands indicating a good photometric agreement between the 
two measurements. 


C.2. WISE: Comparison to other MIR measurements 

WISE photometry can be compared to Spitzer IRAC measure¬ 
ments for the first two bands, to IRAS 12 pm and to MIPS 24 
pm. WISE, IRAC, IRAS and MIPS are not calibrated on the 


same reference spectral shape, so we first apply colour correc¬ 
tions to our flux densities before the comparison. These colour 
corrections are of the order of < 1%, 3%, 4%, 2%, 4%, 5%, <1% 
and 9% for WISE 3.4 pm, WISE 4.6 pm, WISE 12 pm, WISE 
22 pm, IRAC 3.6 pm, IRAC 4.5 pm, IRAS 12 pm, and MIPS 
24 pm respectively. We remove from the comparison all of the 
upper limits in any of the bands and we are left with 35 galaxies 
for WISEl, 37 galaxies for WISE2, 9 galaxies for WISE3 and 
32 galaxies for WISE4. 

The ratios of WISE flux densities to the corresponding MIR 
flux densities are shown in Eig. C.l The median ratios are 
FwISE3.4/FiraC3.6 = 0.97 + 0.12, EwiSE4.6/FiRAC4.5 = 0.87 + 0.17, 
Fwisei2/Firasi2 = 0.88 + 0.13, and Ewise22/Fmips24 = 0.79 + 
0.22. We note that the dispersion is larger for the first two WISE 
bands for WISE flux densities less than ~ 1 mJy. WISE 22 pm 
flux densities are also systematically lower than MIPS 24 pm 
flux densities, reflecting the difference in wavelength between 
the two bands and the steeply rising MIR continuum observed in 
the SED of dwarf galaxies. 

Eor consistency, we also compare the WISE photometry 
from the AllWISE database to our measurements on the maps. 
We confirm that for the brightest and most extended sources, the 
photometry provided by the database underestimates the emis¬ 
sion, and especially at 12 and 22 pm. 


Appendix D: DGS IRS spectra 

D.1. IRS data reduction 


Two versions of the spectra are available in the CASSIS database 
( Lebouteiller et al.|201 1 1 for the staring observations: an “opti¬ 
mal” extraction better suited for point sources and a “tapered 
column” extraction, better suited for partially extended sources. 
A message advises the user if the “tapered column” extraction 
or the “optimal” extraction method should be chosen. The cho¬ 
sen extraction for each of the DGS targets, as well as AORkeys, 
are summarised in Appendix [B| CASSIS uses the AdOpt algo¬ 
rithm (Advanced Optimal extraction, [Lebouteiller et al.]|2010|l 
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Table B.l. DGS SpitzerflRAC and IRS Observing Log. 



IRAC 

IRS 


Sources 

AOR key 

AOR key 

Extraction 

Haro 11 

4326400 

9007104 

Optimal 

Haro 2 

5539840 

9489920 

Map 

Haro 3 

11180288 

12556288 

Tapered 

He 2-10 

4329472 

4340480 

Tapered 

HS 0017+1055 

26387200 

26393344 

Optimal 

HS 0052+2536 

26387456 

17463040 

Optimal 

HS 0822+3542 

4328960 

1763808 

Optimal 

HS 1222+3741 

17564928 

26393600 

Optimal 

HS 1236+3937 

26387712 

26393856 

Optimal 

HS 1304+3529 

26387968 

26394112 

Optimal 

HS 1319+3224 

26388480 

26394624 

Optimal 

HS 1330+3651 

26388736 

26394880 

Optimal 

HS 1442+4250 

10388480 

12562944 

Optimal 

HS 2352+2733“ 

26388992 

26395136 

- 

IZw 18 

4330752 

16205568 

Optimal 

IC 10 

4424960 

26396672 

Map 

II Zw 40 

4327936 

9007616 

Optimal 

Mrk 1089 

11250432 

26395392 

Tapered 

Mrk 1450 

4334336 

16206080 (SL), 
9011712 (LL) 

Optimal 

Mrk 153 

4333056 

4342272 

Optimal 

Mrk 209 

22556672 

12557568 

Optimal 

Mrk 930 

4338944 

4344320 

Tapered 

NGC 1140 

4327168 

4830976 

Tapered 

NGC 1569 

4434944 

3856640 

Tapered 

NGC 1705 

5535744 

9513216 

Map 

NGC 2366 

4436480 

21920768 

Map 

NGC 4214* 

4457984 

- 

- 

NGC 4449 

4467456 

26396928 

Map 

NGC 4861“ 

4337408 

- 

- 

NGC 5253 

4386048 

4386304 

Map 

NGC 625“ 

22520064 

- 

- 

NGC 6822*-“ 

5507072 

- 

- 

Pox 186 

26389248 

12629760 

Optimal 

SBS 0335-052 

4327424 

8986880 

Optimal 

SBS 1159+545 

26389504 

9008896 

Optimal 

SBS 1211+540 

26389760 

26395392 

Optimal 

SBS 1249+493 

26390016 

26395904 

Optimal 

SBS 1415+437 

10392832 

12562432 (SL), 
8990464 (LL) 

Optimal 

SBS 1533+574 

17563904 

8996352 

Optimal 

Tol 0618-402 

4328448 

8090624 

Optimal 

Tol 1214-277 

4336384 

9008128 

Optimal 

UGC 4483 

4329728 

26396160 

Optimal 

UGCA 20“ 

26390272 

26396416 

- 

UM 133 

26390528 

21922304 

Map 

UM 311* 

10392576 

- 

- 

UM448 

4334592 

4342784 

Tapered 

UM461 

4335104 

16204032 (SL), 
9006336 (LL) 

Optimal 

VII Zw 403 

4334080 

9005824 

Tapered 


Notes : 

For these galaxies, the IRS slits are not centred on the source position 
and thus we do not present any IRS spectrum. 

For these galaxies, only local pointings were performed and we 
cannot present an IRS spectrum for the total galaxy. 

Only high resolution IRS spectra (SH and/or LH) are available. 


within the Spectroscopic Modeling Analysis and Reduction Tool 
(SMART, Higdon et al.||2004 i. AdOpt enables optimal extrac¬ 
tion of spectra using a super-sampled PSF. The pipeline includes 
such steps as image cleaning, individual exposure combination 


and background subtraction. Specific attention is given to iden¬ 
tification and removal of bad pixels and outlier rejection at the 
image and spectra levels. 

Some discrepant fluxes between orders at some wavelengths 
have to be removed from the spectra. These cut wavelengths are 
given in Table D.l For the extended galaxies, there is a conver¬ 
sion step from MJy/sr to Jy where the IRS spectrum is multiplied 
by the area over which it has been extracted. The spectra are then 
rescaled to match the observed photometry. 


Table D.l. IRS cut-off wavelengths. 


Module 


Cut-offs (pm) 

SLl 


■ 7.53, 14.02 - A,„,A 

SL2 


- 5.23, 7.49 - 4™,] 

SL3 


LLl 

L'^min 

- 20.5, 37.4 - 4™,] 

LL2 


■ 14.0, 20.52 - A„„A 

LL3 


- 19.8, 21.5-A, J 


D.2. Rescaling the spectra to the photometry 

We describe now in more details the rescaling step of the IRS 
spectrum, and how we derive the correction factor for the SL 
and LL parts of the IRS spectrum. 

The correction factor for LL, Cll can be written; 


ClL - F2aIFirs{M2a) 


(D.l) 


where F 24 is the observed 24 pm MIPS flux density (from Bendo 
et al. ( |2012 i) and Fjrs (M 24 ) is the synthetic photometry for the 
IRS LL, spectrum at 24 pm. Eq. D.l is adapted for WISE 22 pm 
when MIPS 24 pm is not available. 

Several constraints are available for the SL module, and 
used simultaneously to derive a correction that smoothly de¬ 
pends on wavelength, Csl('I)- We fit a spline function, using the 
IDL SPLINE_P procedure, to Csl(I 3), Csi(I4), Csl(W3) and Cn. 
We use Cll here to have a better constraint at the end of the 
spline. We impose a derivative of 0 at 24 pm in order to have 
a smooth function. We also assume that the correction factor 
Csl(3) should always be > 1.0, as the photometry measurement 
always encompasses more flux than the IRS slit. 

Csl(I3), Csl(I4), Csl(W3) are given by: 


( Csl(I^) - F^^^/FjRsih.s) 

1 Csl(I 4 ) = Fs.o/F,Rs(h.o) (D.2) 

( Csl(W 3 ) - F12 IF IRS,sl(Wi 2 ) 

where F^^, F^ q and F 12 are the observed 5.8 and 8.0 pm IRAC 

and the 12 pm WISE flux densities (from Tables]^ and 0; and 

FiRsIh.s), FiRs(h.o) and Fjrs.slIWu) are the synthetic photom¬ 
etry for the IRS SL spectrum at 5.8, 8.0 and 12 pm. 

Note that there is an overlap between the WISE 12 pm Al¬ 
ter and the LL wavelengths. As the LL spectrum has already 
been corrected, we must And the flux missing in the SL spec¬ 
trum, i.e., Firs,sl(Wi 2 ), to match F12. But the integration over 
the WISE bandpass is not linear, i.e.,: Firstsl(Wi 2 ) + Firs(W 12 ) 
- F IRS,ll( 331 u)■ Thus there is no simple way of deriving Csl(3313). 
Instead, we apply the following method; 


1. We generate a grid of potential Csl(W3), (x), from 0.01 to 

10 ; 
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Fig. D.l. Example of IRS treatment for HS 1304+3529; the raw 
spectrum is presented in the top panel. The photometry points 
are overlaid in filled purple diamonds as well as synthetic pho¬ 
tometry as blue open diamonds. The galaxy is not detected at 
5.8 pm as indicated by the open purple triangle. Green, orange 
and red parts of the spectrum represents the SL2 and LL2, SL3 
and LL3, SLl and LLl spectra, respectively. The middle panel 
shows the spectra after cutting-off the discrepant wavelengths 
and rescaling to match the photometry. The spline used to rescale 
SL is shown in Fig. |D.2| The bottom panel shows the spectra 
smoothed to a S/N of 3. 


2. find the spline going through Csl(I3), Csl(I4), Cll, and each 
X/, sp(x;); 

3. correct the SL spectrum with each sp(x,) and 

4. compute the synthetic WISE 12 pm photometry for each cor¬ 
rected total IRS spectrum, FiRs(W\ 2 )(sp(xi)). 


HS1304+3529 



Wavelength [/xm] 

Fig. D.l. Example of a spline used to rescale SL in the case of 
HS 1304+3529. The diamonds represent the ratio between the 
observed and the synthetic photometry. The spline has been fit¬ 
ted to the blue diamonds, except at 5.8 pm where the galaxy is 
not detected. The black curve is the final adopted spline, sp(xo), 
where ;ro = L9. The orange crosses are the corrections applied 
to the SL part of the spectrum. 


In the case of NGC 1140 and NGC 1569, the IRS SL slit 
only covers a small part of these extended galaxies. We use the 
ISOCAM spectrum (from Galliano et al.||20()3 2005| l for com¬ 
parison after applying the same rescaling step. The two spec¬ 
tra are consistent with each other and with the general NIR - 
MIR shape of the SED for both galaxies. Thus we consider that 
these IRS spectra are reliable and use them to constrain the SED 
model. 


Appendix E: Specific SED modeiiing 

In this Appendix, we detail the specific SED modelling for sev¬ 
eral DGS sources, and show the SEDs for the galaxies with and 
without any additional features (addition of a MIR modBB, in¬ 
cluding extinction, etc.) to show how the fit is improved. 

- Haro 11 - As this galaxy has evident silicate absorption fea¬ 
tures in the IRS spectrum around 10 pm, we add extinction in 
our modelling. The output luminosity is modified by: 

Lf ■‘'“'(T, Ay) = Lf •“'“'(d. Ay = 0) X e-^oWA./1.086 (g 


Csl(W3) is the xq that gives FjRs(Wi 2 )(sp(xo)) = Fn. 


The final IRS spectrum, Firs,c 
given by: 


r (Asl) and Firs, corr(ALL), is 


FlRSyorriAsL.) - Fjrs (As l) X Cs l(Asl) 
FIRSyorriAlL) - F jrs (Ah) X Cll 


(D.3) 


The method presented here is adapted depending on the 
number of constraints for each galaxy. Upper limits are not con¬ 
sidered for the correction of IRS spectrum. An example of this 
treatment of IRS spectrum is shown in Fig. D.l along with the 
spline used to correct the SL spectrum in Fig. D.2 The final IRS 
spectra for the DGS galaxies are shown in Fig.|D.3| 


where tq is the optical depth, determined from the optical 
properties and the size distribution of our dust composition. The 
extinction parameter Ay is left varying between 1 and 20, with a 
best fit value of Ay = 5. The quality of the fit is clearly improved 
in the NIR-MIR range after introducing Ay (see Fig. ED- 
Because the stellar continuum is also affected by extinction, the 
stellar template used for the fit does not provide a reliable fit of 
the 2MASS NIR data anymore. Instead, we replace the stellar 
component by a blackbody and impose a temperature T > 1000 
K. We obtain a best fit temperature T- 30,000 K, corresponding 
to a maximum of the blackbody emission around 0.1 pm in the 
UV - a sign of a young stellar population. This is consistent 
with the results of |Adamo et al.| ( |2010| l who found a maximum 
in the cluster formation history ~ 3.5 Myr ago, and with the UV 
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Fig.D.3. IRS spectra for the DGS galaxies. The SL module is shown in green and the LL module is shown in red. The positions of 
the [SIV] 10.5 pm, [Nell] 12.8 pm, [Nelll] 15.6 pm, [Sill] 18.7 pm and 33.5 pm spectral lines are indicated. Line intensities have 
been extracted by Cormier et al. (|2015|l. 
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Fig. D.3. IRS spectra for the DGS (continued). 
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Fig. D.3. IRS spectra for the DGS (continued). 


part of the SED presented in Cormier et al. (20121. 


- HS 0822+3542 - The observed PACS upper limits at 70 and 
100 pm are not consistent with the rest of the MIR to FIR 
photometry (see Fig. El- We recommend using the synthetic 
photometry provided by the model as upper limits for the 70 
and 100 pm PACS wavelengths, i.e., F^o < 41 mJy and Fmo < 
48 mJy. These values are reported in Table 


- HS 2352+2733 - There are only 5 detected points for this 
galaxy, and thus not enough constraints to ht a full SED. 
However, we can get a rough estimate of the dust mass by htting 
a modBB with a fixed (3-2.0, through MIPS 24 pm, PACS 70 
pm and PACS 100 pm where the galaxy is detected. We obtain 
a temperature of T= 52 K and M^b = 1.01 x 10^ Mq. This mass 
can be seen as a lower limit to the real dust mass (see Sect. |3.4| l. 
This value is reported in Table and marked with a 


- II Zw 40 - This galaxy has an IRS spectrum with very good 
S/N. However, when compared to other MIR measurements, 
it seems that the Spitzer IRS-LL spectrum is not consistent 
with the rest of the MIR - FIR photometry, meaning that for 
this galaxy, a wavelength dependent correction may also be 
needed for the IRS LL spectrum. As an alternative, we use 


the ISOCAM spectrum for this galaxy (from Galliano et al. 


2005|l, which covers a 5.6 to 16.3 pm wavelength range. The 


ISOCAM spectrum was matched to the broad band photometry 
by applying the same rescaling process as for the IRS SL 
spectrum. To properly match the end of the ISOCAM spectrum, 
we add an extra modBB i(3 fixed to 2.0 and best fit temperature 
113 K). 


- Mrk 153 - This galaxy has a prominent silicate emission 
feature in its IRS spectrum around 10 pm. This emission 
feature originates from hot small silicate grains, possibly in the 
accretion disk around an AGN. We allow the the silicate-to- 


(silicate+graphite) grain mass fraction to vary in our ht and get 
a ratio of 0.94 (i.e., 1.3 times the Galactic value). 


- SBS 0335-052 - This galaxy has a very surprising SED with 
a IR peak around 15 - 30 pm, and an IRS spectrum showing 
silicate absorption features superimposed on a featureless 
continuum (see also Thuan et al.|1999 [Dale et al.|2001t |Houck| 


et al.|2004| Galliano et al.|2008 1 . The silicate absorption feature 

around 10 pm observed in the IRS spectrum indicates that 
extinction is present in this galaxy. We add extinction in our 
model as for Haro 11, and get Ay =1.1. The peculiar shape of 
the IRS spectrum requires two MIR modihed blackbodies to 
obtain a satisfactory ht (J3 hxed to 2.0 and best ht temperatures 
Ti = 273 K and T 2 = 123 K, see Fig.[ET]i. 


- Tol 1214-277 - The MIPS 70 pm point is not consistent with 
the rest of the photometry; IRS and MIPS 24 pm on one side 
and PACS 70 and 100 pm on the other side (see Fig. [T 2 I 1 . It is 
another conhrmation that this point is discrepant (|Remy-Ruyer 


et al.|2013 1, and thus we do not consider it in the modelling. 
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Fig.E.l. Final SEDs in black compared to the SED without the additional features in blue (MIR modBB, extinction, see text 
source per source details). Note how the final SEDs better reproduce the observed IR emission for these galaxies. 





















